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Abstract 

An  imaging  technique  called  indirect  photography  was  recently  demonstrated  for 
imaging  otherwise  hidden  scene  information  through  the  collection  and  radiometric 
modeling  of  light  reflecting  off  of  visible  reflecting  surfaces.  A  laser  is  reflected  from 
a  visible  surface  onto  the  hidden  object.  A  camera  then  collects  light  which  reflects 
from  the  object  back  to  the  visible  surface.  The  camera  images  are  used  to  construct 
an  indirect  image  of  the  hidden  object.  The  ability  to  focus  the  beam  that  reflects 
off  the  visible  surface  to  a  tight  spot  on  the  hidden  object  is  theorized  to  simulate 
a  direct  line  of  sight  from  the  light  source  to  the  object,  effectively  reducing  the 
radiometric  model  to  that  of  the  previously  demonstrated  dual  photography,  and 
leading  to  much-simplified  results.  This  work  demonstrates  that  capability. 

Experiments  involving  the  focusing  of  light  transmitted  through  highly  scatter¬ 
ing  media  have  been  adapted  to  the  case  of  reflection.  The  focusing  of  a  HeNe 
laser  reflected  off  diffusely  scattering  surfaces  to  a  tight  spot  is  demonstrated  using 
phase  control.  A  liquid  crystal  spatial  light  modulator  assigns  phase  delay  pixel  by 
pixel  to  the  Gaussian  beam  to  conjugate  the  scattering  properties  of  the  reflector. 
Intensity  enhancement  values  between  13.8  and  122.3  are  achieved  for  the  surfaces 
studied.  While  samples  with  relatively  strong  specular  reflection  achieve  higher  en¬ 
hancement  for  a  given  geometry,  the  more  diffusely  reflecting  materials  achieve  higher 
enhancement  in  regions  of  nonspecular  scatter.  High  diffuse  scatter  thus  facilitates 
the  non- mechanical  steering  of  the  focused  spot. 
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THE  FOCUSING  OF  LIGHT  SCATTERED  FROM  DIFFUSE  REFLECTORS 


USING  PHASE  MODULATION 

I.  Introduction 


1.1  Research  Motivation 

The  aim  of  this  thesis  is  to  focus  light  scattered  off  of  diffuse  surfaces  by  adjusting 
the  phase  elements  of  the  beam  with  a  liquid  crystal  (LC)  spatial  light  modulator 
(SLM).  This  focusing  capability  is  theorized  to  have  implications  for  advancing  ef¬ 
forts  in  an  imaging  technique  called  “indirect  photography”  The  concept  of  indirect 
photography  has  recently  been  introduced  as  a  photographic  technique  for  imaging 
objects  hidden  around  corners.  Indirect  photography  stems  from  a  technique  called 
“dual  photography,”  a  method  which  exploits  Helmholtz  reciprocity  by  using  a  light 
source  to  illuminate  a  scene  and  measuring  the  light  transport  between  each  pair 
of  source-to-camera  pixels.  The  “dual  image”  is  generated  with  the  computational 
interchanging  of  the  source  and  camera.  Attaining  an  “indirect  image”  involves  col¬ 
locating  the  camera  and  the  light  source  to  recover  parts  of  the  scene  not  directly 
visible  to  either  the  camera  or  the  controlling  light  source  [1] . 

Collocating  the  light  source  and  the  camera  gives  rise  to  a  potentially  powerful 
imaging  technique  which  could  become  a  valuable  asset  in  U.S.  military  operations.  In 
effect,  a  soldier  could  “see  around  corners”  by  imaging  an  object  with  light  reflected  off 
of  a  diffuse  surface.  This  technique  uses  the  radiometric  principles  developed  in  dual 
photography  while  adding  another  diffuse  surface  reflection  to  factor  into  the  imaging 
process  [2],  The  ability  to  focus  the  light  source  onto  the  object  after  reflection  off  of 
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the  first  diffuse  surface  is  assumed  to  be  equivalent  to  using  a  light  source  with  a  direct 
line  of  sight  to  the  object,  thereby  reducing  the  problem  of  indirect  photography  to 
dual  photography. 

Recent  experiments  by  Vellekoop  et  al.  have  shown  that  on  transmission,  solving 
for  and  applying  a  uniquely  shaped  wavefront  can  be  used  to  steer  light  through 
a  highly  scattering  object  in  a  process  called  “inverse  diffusion.”  The  wavefront  is 
shaped  to  exactly  match  the  scattering  properties  of  the  object  so  that  the  object 
focuses  the  light  to  a  point  [3] .  Inspired  by  this  research,  the  method  introduced  here 
to  focus  light  scattered  upon  diffuse  reflection  parallels  the  technique  developed  for 
inverse  diffusion  in  transmission.  For  the  purpose  of  this  thesis,  we  will  refer  to  inverse 
diffusion  as  demonstrated  by  Vellekoop  as  “transmissive  inverse  diffusion”  while  using 
the  term  “reflective  inverse  diffusion”  to  signify  the  adaptation  of  the  technique  to 
the  case  of  reflection. 

1.2  Problem  Statement 

This  thesis  seeks  to  accomplish  a  set  of  four  overarching  goals  including  (1)  the 
demonstration  of  a  proof  of  concept  for  focusing  light  scattered  off  of  an  assortment  of 
diffusely  reflecting  surfaces  with  various  scattering  properties,  (2)  the  nonmechanical 
steering  of  the  focused  spot,  (3)  the  identification  of  certain  material  and  experimental 
factors  which  must  be  considered  when  applying  reflective  inverse  diffusion  to  indirect 
photography,  and  (4)  a  prescriptive  outline  sketching  out  the  progression  required 
to  use  the  results  from  this  thesis  and  extend  them  to  practical  implementation  in 
indirect  photography. 
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1.3  Methodology 


The  optical  setup  used  for  demonstrating  reflective  inverse  diffusion  is  designed 
to  contain  aspects  of  Vellekoop’s  experiment  for  transmissive  inverse  diffusion  while 
removing  constraints  which  would  render  it  unuseful  for  indirect  photography.  Optical 
components  such  as  a  microscope  objective  placed  shortly  before  and  shortly  after  the 
scattering  sample  are  removed.  The  technique  is  demonstrated  on  six  nonspecular 
reflective  surfaces  and  the  differences  in  the  performance  of  the  technique  are  analyzed 
with  respect  to  material  properties.  The  main  idea  of  the  experiment  is  to  achieve 
focused  spots  from  diffuse  reflection  off  of  various  reflecting  materials  while  using  a 
feasible  indirect  photography  setup,  and  to  steer  the  focused  spot  to  different  locations 
on  a  simulated  object  by  choosing  different  target  areas  for  the  same  setup. 

1.4  Overview 

This  thesis  is  broken  up  into  five  chapters.  The  second  chapter  discusses  back¬ 
ground  from  literature  on  dual  photography,  indirect  photography,  radiometry,  and 
the  properties  of  the  surfaces  used  in  this  experiment.  Chapter  3  describes  the  test¬ 
ing  methodology.  Chapter  4  details  the  results  of  the  experiments,  and  Chapter  5 
offers  recommendations  for  improving  the  technique  used  in  this  thesis  and  for  further 
developing  it. 
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II.  Literature  Review 


2.1  Introduction 

This  chapter  reviews  previous  research  which  shapes  and  applies  to  the  technique 
developed  in  this  thesis  for  focusing  light  scattered  from  a  diffuse  surface.  While  this 
technique  does  not  require  detailed  knowledge  of  the  theoretical  underpinnings  of 
dual  and  indirect  photography,  a  brief  review  of  the  theory  and  literature  are  offered 
to  contextualize  the  problem  and  experimental  setup.  The  principles  of  transmis¬ 
sive  inverse  diffusion  on  which  the  reflective  inverse  diffusion  technique  is  based  are 
discussed.  Nonmechanical  beam  steering  using  a  LC  SLM  is  also  introduced  and 
related  to  beam  steering  via  reflection  off  of  highly  scattering  surfaces.  Finally,  char¬ 
acteristics  of  the  diffuse  surfaces  chosen  to  demonstrate  reflective  inverse  diffusion  are 
highlighted. 

2.2  Dual  Photography 

Dual  photography  as  developed  by  Sen  et  al.  is  a  completely  image-based  math¬ 
ematical  method  for  generating  pictures  from  the  viewpoint  of  a  light-emitting  pro¬ 
jector  rather  than  a  camera.  A  pixelated  light  source  is  used  to  illuminate  a  scene 
while  a  digital  camera  is  used  to  record  images  from  reflection  of  the  scene  or  from 
reflection  off  of  a  diffuse  surface.  The  images  are  used  to  construct  a  transport  ma¬ 
trix  T  which  maps  light  source  pixel  to  camera  pixel  transport  characteristics.  As 
illustrated  on  the  left  in  Figure  1,  the  linearity  of  light  transport  allows  for  a  matrix 
mapping  between  source  and  camera  given  by 

c'  =  Tp'  (1) 
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in  which  p'  is  an  mnxl  vector  representing  the  light  source  with  mxn  pixels,  d  is  a 
pqxl  vector  representing  a  camera  with  pxq  pixels,  and  T  is  the  transport  matrix  of 
size  pqxmn  which  tracks  each  possible  transport  path  from  source  pixel  to  camera 
pixel  [2], 


Figure  1.  In  the  primal  configuration,  light  is  emitted  by  a  real  projector  and  captured 
with  a  real  camera.  The  matrix  T  describes  the  light  transport  from  projector  to 
camera,  where  T(J  is  the  transport  coefficient  from  projector  pixel  j  to  camera  pixel 
i.  In  the  dual  configuration,  the  camera  and  projector  are  reversed,  and  Helmholtz 
reciprocity  specifies  that  Y  j  ■  =  /  .  meaning  T"  =  Tt  . 


As  the  right  side  of  Figure  1  indicates,  switching  the  positions  of  the  source  and 
camera  is  achieved  mathematically  by  recognizing  the  light  transport  characteristics 
between  each  set  of  pixels  remain  unchanged,  and  the  transpose  of  the  transport 
matrix  can  be  used  to  solve  for  an  image  of  the  scene  from  the  perspective  of  the 
projector  in  accordance  with 

p  =  Ttc".  (2) 

The  user  can  manipulate  images  of  the  scene  with  a  variety  of  simulated  light  source 
characteristics  making  this  a  particularly  useful  technique  in  attaining  photorealism 
in  computer  graphics  applications.  Sen  et  al.  performed  an  experiment  using  dual 
photography  to  image  the  illuminated  face  of  a  playing  card  upon  its  reflection  off  the 
page  of  a  book.  While  the  camera  did  not  have  access  to  the  face  of  the  card,  the  card 
was  set  so  that  the  only  light  that  reached  the  camera  had  undergone  diffuse  bounces 
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at  the  card  and  at  the  book.  The  transport  matrix  was  acquired  from  5742  images  to 
produce  a  dual  image  with  resolution  66x87  from  the  perspective  of  the  projector  [4], 
Taking  away  the  constraint  that  the  projector  or  light  source  requires  a  direct  line  of 
sight  to  the  object  being  imaged-the  playing  card,  in  this  case-would  give  rise  to  an 
imaging  technique  which  would  serve  as  a  powerful  capability  in  military  operations. 
Figure  2  shows  the  relationship  between  a  dual  setup  and  an  indirect  setup. 


Figure  2.  Real  world  application  of  (a)  dual  photography  and  (b)  indirect  photography. 

2.3  Indirect  Photography 

Indirect  photography  is  an  extension  of  dual  photography  which  removes  the  re¬ 
quirement  that  the  light  source  has  a  direct  line  of  sight  to  the  object  being  imaged. 
The  technique  involves  collocating  the  camera  and  a  HeNe  light  source.  The  most 
recent  set  of  experiments  used  the  setup  shown  in  Figure  3,  for  which  the  camera  and 
laser  are  placed  without  a  direct  line  of  sight  to  the  checkerboard  underneath  the  top 
surface.  The  laser  hits  a  spot  on  the  first  surface  (aq,  y\)  and  is  reflected  onto  the  top 
surface  (x2,  y2,  z2)  where  the  checkerboard  sits.  Light  which  is  reflected  off  the  top 
surface  and  onto  the  first  surface  within  the  camera’s  field  of  view  (.1:3,  y3)  is  collected 
by  the  camera. 

For  both  dual  and  indirect  photography,  the  final  reflecting  surface  is  imaged.  The 
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Figure  3.  Most  recent  experimental  setup  for  indirect  photography.  The  laser  reflects 
off  the  first  surface  at  the  points  labeled  (x\,yi),  the  object  to  be  imaged  which  is  here 
a  2x2  checkerboard  with  coordinates  labeled  (^2  >2/2  5-22  )>  and  finally  back  to  the  first 
surface  at  the  points  labeled  (#3 ,2/3).  The  first  and  last  surface  are  coplanar  with  the 
latter  being  defined  by  the  camera  field  of  view. 


application  of  radiometric  theory  and  a  blind  deconvolution  process  to  these  images 
yields  an  image  of  the  object.  For  a  thorough  detailing  of  the  radiometric  development 
involved  in  attaining  dual  and  indirect  images,  see  [2]. 

2.4  Inverse  Diffusion  using  an  SLM 

The  ultimate  goal  of  this  thesis  is  to  provide  a  framework  for  simplifying  and 
improving  the  indirect  photography  technique  through  focusing  light  to  a  point  after 
reflection  off  a  diffuse  surface.  The  ability  to  focus  to  a  point  on  reflection  effectively 
reduces  the  problem  of  indirect  photography  to  that  of  dual  photography  by  removing 
the  complexity  which  comes  with  reflecting  the  light  source  off  the  first  diffuse  surface. 
We  assume  that  focusing  upon  reflection  off  of  a  diffuse  surface  to  a  point  on  the 
object  is  equivalent  to  illuminating  the  object  with  a  source  with  a  direct  line  of 
sight.  We  conjecture  that  focusing  upon  reflection  can  be  done  in  a  manner  similar 
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to  the  focusing  of  light  in  transmission  by  way  of  inverse  diffusion  using  an  SLM  and 
that  many  of  the  same  governing  relationships  apply. 

Transmissive  inverse  diffusion  is  a  technique  by  which  an  LC  SLM  is  used  to 
shape  a  wavefront  to  focus  light  through  opaque  objects.  When  an  unshaped  beam 
is  transmitted  through  a  disordered  scattering  medium,  the  result  is  a  complicated, 
random  interference  pattern  of  many  waves  forming  laser  speckle.  Using  the  LC  SLM, 
the  beam  is  broken  into  hundreds  or  thousands  of  light  beams  with  phase  delays  with 
respect  to  the  each  other.  The  resulting  total  field  at  a  point  behind  the  sample  is 
the  sum  of  the  speckle  patterns  of  these  individual  beams  [3].  Figure  4  illustrates 
how  shaping  the  beam  which  impinges  on  the  scattering  sample  serves  to  alter  the 
’’random  walk”  of  the  light  through  the  sample.  The  end  result  is  the  constructive 
interference  of  transmitted  light  at  the  designated  focal  point. 


Figure  4.  Opaque  lensing.  a)  A  plane  wave  impinging  on  a  scattering  object  performs 
a  random  walk  through  the  material.  The  light  that  makes  it  through  is  scattered  in 
all  directions,  b)  The  wave  is  shaped  to  match  the  scattering  in  the  material  so  that 
the  object  focuses  the  wave  to  any  designated  point. 


Assigning  phase  delay  to  induce  speckle  patterns  which  align  to  one  focused  spot, 
the  total  field  Em  in  the  target  area  is  given  by 


where  Nt  represents  the  total  number  of  modulated  segments  of  the  SLM,  trnn  is  the 
transmission  coefficient,  and  An  and  <j)n  are  the  amplitude  and  phase  of  the  laser  light 
reflected  from  the  SLM  segment  n.  The  constants  tmn  are  statistically  independent 
and  follow  a  circular  Gaussian  distribution,  leading  to  an  expression  for  the  average 
amplitude  enhancement  (a)  given  by 


<«>  = 


(  E0pt\)r 


7 r 


—  A/  4  (nt  -  1)  +  1 


\Nt- 


(4) 


( |  Ernd  | ) 

where  (\Eopt\)rms  is  the  RMS  value  for  the  field  amplitude  within  the  optimized  target 
area  and  ( |  Ern(i  | )  rms  is  the  RMS  field  amplitude  of  the  random  speckle  pattern  before 
optimization  begins.  The  maximum  intensity  enhancement  rj  =  a2  is  the  ratio  be¬ 
tween  the  optimized  intensity  and  the  average  intensity  of  the  speckle  pattern  before 
optimization  [5]. 

Comparing  the  final  optimized  intensity  to  the  surrounding  speckle  intensity  is  of 
greater  consequence  when  considering  its  role  in  indirect  photography.  The  goal  is  to 
focus  to  a  point  on  the  object  and  minimize  contributions  from  light  reflecting  off  of 
the  object  outside  the  focus.  The  reflective  inverse  diffusion  enhancement  is  thus 
defined  as  the  ratio  between  the  optimized  intensity  and  the  average  intensity  of  the 
speckle  pattern  surrounding  the  focus. 

To  account  for  the  time-dependent  decorrelation  of  the  speckle  pattern  due  to 
environmentally-induced  instability,  Vellekoop  estimates  that  the  effective  enhance¬ 
ment  decreases  according  to 


V  ^ 

Veff  -i  ,  NtT  (5) 

^  Ts 

where  T  is  the  time  required  for  a  single  measurement  and  Ts  is  the  time  scale  for 
which  the  speckle  pattern  of  the  sample  remains  stable.  The  persistence  time  is 
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defined  as  the  decay  time  of  the  transmitted  speckle  autocorrelate  and  depends  on 
both  the  type  of  sample  used  and  on  on  the  environmental  conditions.  The  persistence 
time  of  the  paint  sample  Vellekoop  presents  is  reported  to  be  90  minutes  due  primarily 
to  humidity  in  the  sample  [6].  The  persistence  time  of  the  speckle  patterns  from 
reflection  off  of  diffuse  surfaces  is  expected  to  be  higher  than  for  transmission  through 
highly  scattering  materials.  In  the  latter,  the  beam  undergoes  hundreds  of  scattering 
events  within  the  sample,  requiring  that  the  sample  is  stable  with  sub- wavelength 
accuracy  throughout  the  optimization  process  [3]. 

2.5  Liquid  Crystal  Spatial  Light  Modulator 

The  LC  SLM  is  becoming  a  highly  researched  and  tested  optical  component  for 
military  imaging  sensors  as  a  result  of  its  small  size,  light  weight,  and  adaptability  [7], 
making  it  a  good  candidate  for  its  potential  role  in  indirect  photography.  The  LC 
devices  used  in  these  experiments  are  of  the  Nematic  variety,  which  are  uniaxially  bire- 
fringent.  Applying  an  external  AC  electric  field  induces  a  dipole  in  the  LC  molecule 
to  produce  a  torque  and  resulting  rotation  of  the  optical  axis.  The  incident  wave- 
front  experiences  an  optical  path  difference  (OPD)  based  on  the  orientation  of  the  LC 
molecule.  The  OPD  can  be  adjusted  on  a  pixel-by-pixel  basis,  typically  achieving  up 
to  a  2tt  phase  depth.  Wavefront  control  is  thus  accomplished  modulo-27r[8].  Figure  5 
shows  a  top- view  schematic  of  a  reflective  LC  SLM. 

2.6  Nonmechanical  Beam  Steering 

To  further  aid  in  simplifying  indirect  photography,  the  ability  to  nonmechanically 
steer  the  focused  spot  would  allow  for  scanning  the  hidden  object  to  be  imaged 
without  adjusting  the  position  of  any  optical  elements.  When  working  with  a  coherent 
beam  without  an  intermediate  diffuse  reflection,  applying  a  saw-tooth  phase  profile 
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VLSI  Die 

PIN  GRID  ARRAY  PACKAGE 


Figure  5.  Top  view  schematic  of  SLM.  Light  is  reflected  by  the  electrode  mirrors  of 
the  pixels  and  experiences  phase  delay  in  the  liquid  crystal.  The  VLSI  die  and  pin  grid 
array  package  control  the  voltage  differences  which  align  the  liquid  crystal  molecules 
according  to  the  desired  index  of  refraction  and  resulting  phase  delay  at  each  pixel. 


to  the  SLM  effectively  induces  an  effect  similar  to  a  blazed  diffraction  grating.  The 
wavefront  is  tilted  and  deflected  to  a  nonzero  angle  as  illustrated  in  Figure  6.  A 
more  significant  steer  angle  is  achieved  with  narrower  phase  ramps  ranging  from  0 
to  2ir  phase  delay.  Since  the  SLM  has  a  finite  number  of  addressable  pixels  and 
the  smooth  ramp  becomes  more  noticeably  discretized  with  fewer  intermediate  phase 
delay  values,  the  steer  angle  is  limited. 


Figure  6.  Beam  steering  schematic.  The  tilted  wavefront  (dashed  red)  results  from 
application  of  a  modulo  2n  saw-tooth  phase  pattern.  A  greater  tilt  is  induced  with 
narrower  modulated  columns,  which  is  limited  by  the  number  of  addressable  pixels  of 
the  SLM. 


The  steer  angle  9  of  a  beam  using  this  technique  is  geometrically  described  with 
6  =  arcsin  ^  where  7  represents  the  pitch  of  the  phase  ramp  and  A  is  the  wavelength 
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[9].  The  displacement  d  of  the  beam  is  therefore  limited  geometrically  by  the  focal 
length  of  the  lens  used  to  focus  the  beam  onto  the  target  area  according  to  the  relation 
d  =  f  tan  0  m  fO  for  small  values  of  6. 

When  using  the  SLM  to  “steer”  a  focused  spot  in  the  context  of  reflective  inverse 
diffusion,  the  only  option  currently  available  is  selecting  a  new  target  area  for  focus 
and  running  through  the  optimization  process  as  described  in  Section  3.6  for  the  new 
location.  The  limit  to  which  the  focused  spot  is  steered  to  a  given  point  is  a  function 
of  the  intensity  of  scattered  light  at  the  corresponding  scatter  angle.  The  fundamental 
geometric  descriptor  of  reflectance  is  given  by  the  bidirectional  reflectance  distribution 
function  (BRDF)  fr  (0,  .  dp  6b,  ds),  describing  the  reflectance  in  a  certain  direction  as 
specified  by  the  angles  defined  in  Figure  7.  It  is  expected  that  the  ability  to  “beam 
steer”  is  highly  dependent  on  the  BRDF  of  each  particular  sample.  All  steering 
experiments  performed  in  this  thesis  are  described  in  terms  of  these  angles. 


Figure  7.  Schematic  illustration  of  the  bidirectional  reflectance  distribution  function 
(BRDF). 


The  polar  angle  9  is  measured  from  the  surface  normal  while  the  azimuth  angle 
cf)  is  measured  from  an  arbitrary  reference  in  the  surface  plane  [10].  For  the  current 
purpose,  this  reference  is  the  plane  containing  the  incident  beam  (d*  =  0). 
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2.7  Reflective  Inverse  Diffusion  within  Indirect  Photography 

To  tie  together  the  concepts  which  have  been  laid  out,  the  progression  for  which 
reflective  inverse  diffusion  will  theoretically  reduce  indirect  photography  to  dual  pho¬ 
tography  while  imaging  components  still  have  no  line  of  sight  to  the  object  should  be 
outlined.  This  thesis  aims  to  provide  a  proof  of  concept  for  reflective  inverse  diffusion, 
demonstrate  nonmechanical  beam  steering  within  the  context  provided  in  Section  3.9, 
and  establish  a  framework  by  which  enhancement  can  be  predicted  with  knowledge 
of  the  scattering  properties  of  various  materials.  The  arrayed  detector  (AD)  used  to 
detect  the  target  intensity  currently  takes  the  place  of  the  hidden  object. 

The  motivation  to  collocate  the  light  source  with  the  camera  without  a  line  of 
sight  to  the  object  is  analogous  to  the  motivation  to  provide  feedback  to  the  SLM 
from  a  detector  which  is  collocated  with  the  camera:  the  potential  value  of  indirect 
photographing  techniques  lie  with  the  ability  to  image  an  object  which  is  hidden 
from  all  optical  components  which  are  used  to  image  it.  Figure  8  illustrates  how 
introducing  a  second  AD  to  record  intensity  pattern  data  serves  first  to  establish  a 
relationship  between  patterns  of  scatter  at  the  target  and  at  a  location  displaced  from 
the  target;  the  second  AD  subsequently  uses  the  relationship  to  provide  feedback  to 
the  SLM.  Furthermore,  this  is  done  repeatedly  until  the  AD  which  is  gathering  scatter 
information  from  the  diffuse  surface  is  collocated  with  the  laser  and  the  remainder  of 
the  imaging  system. 

2.8  Diffusely  Reflecting  Materials 

Six  reflectors  are  chosen  based  on  the  extreme  variation  of  their  scattering  proper¬ 
ties.  The  differences  in  the  materials’  bidirectional  reflectance  distribution  functions, 
surface  roughnesses,  and  their  resulting  achieved  enhancement  offer  insight  into  po- 


13 


Figure  8.  SLM  feed  back  progression.  To  make  reflective  inverse  diffusion  useful  in 
the  context  of  indirect  photography:  (a)  Gather  intensity  pattern  information  from 
another  arrayed  detector  (AD)  displaced  from  the  focused  spot  and  establish  a  rela¬ 
tionship.  (b)  Establish  a  relationship  between  the  intensity  patterns  and  optimize  using 
feedback  from  the  second  AD.  (c)  Repeat  for  greater  displacements  until  (d)  informa¬ 
tion  gathered  from  an  AD  collocated  with  the  laser  can  be  used  to  control  the  SLM  to 
achieve  a  focused  spot  on  the  first  detector. 


tential  considerations  in  using  this  technique  for  indirect  photography.  The  materials 
used  include  Spectralon,  Infragold,  a  flat  white  paint  sample,  graphite,  and  two  alu¬ 
minum  samples. 

Spectralon  is  a  translucent  material  consisting  of  sintered  and  compressed  polyte- 
trafluoroethylene.  Its  high  reflectivity,  low  absorption,  and  nearly  Lambertian  radia¬ 
tive  properties  make  it  a  useful  radiometric  calibration  material  in  the  infrared  and 
visible  spectra  [11].  Infragold  coating  is  another  highly  Lambertian  material  more 
commonly  used  in  the  infrared  region  of  the  spectrum  for  its  excellent  reflectance 
properties,  but  is  also  relatively  reflective  in  the  632.8  nm  region. 

The  white  paint  sample  consists  of  a  0.25  mm  silica  cover  slip  with  a  thin  layer 
of  white  spray  paint,  similar  to  the  strongly  scattering  sample  used  in  transmissive 
inverse  diffusion.  It  is  used  to  verify  that  the  algorithm  and  equipment  are  set  up 
correctly  in  a  reproduction  of  Vellekoop’s  experiment  [3].  It  is  also  used  on  reflection 
as  another  nearly  Lambertian,  highly  scattering  surface;  however,  with  a  layer  of  paint 
having  a  thickness  on  the  order  of  10  ±4 /rm,  losses  due  to  transmission  decrease  initial 
starting  speckle  intensity. 

Graphite  is  a  highly  absorbing  material  in  the  visible  region  [12],  and  the  two 
aluminum  samples  have  significantly  different  surface  roughnesses.  One  is  lightly 
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brushed  while  the  other  is  sandblasted,  yielding  two  aluminum  samples  with  differing 
BDRF’s.  The  brushed  sample  has  a  stronger  peak  reflectance  in  the  specular  region 
as  evidenced  by  its  shine.  Section  3.8  displays  photos  of  all  the  samples. 

2.9  Summary 

Reflective  inverse  diffusion  is  a  novel  concept  which  is  adapted  from  work  done 
by  Vellekoop  et  al.  on  transmission.  Its  potential  use  in  reducing  the  complexity 
of  indirect  imaging  to  the  imaging  techniques  applied  in  dual  photography  provides 
both  the  research  motivation  and  the  framework  for  determining  a  relevant  geometric 
setup  for  experiments.  The  literature  reviewed  in  this  chapter  provides  background 
and  context  for  the  methodology  presented  in  Chapter  3. 
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III.  Methodology 


3.1  Introduction 

Using  the  background  developed  in  Chapter  2,  the  methodology  for  this  thesis 
is  aimed  towards  providing  a  proof  of  concept  for  reflective  inverse  diffusion  while 
operating  in  framework  prescribed  by  the  goal  of  its  use  in  indirect  photography.  The 
optical  setup  outlined  in  Section  3.4  is  used  consistently  for  a  set  of  diffusely  reflecting 
materials  with  various  scattering  properties. 

3.2  Laboratory  Equipment 

The  main  components  of  the  laboratory  setup  include  a  632.8  nm  HeNe  laser, 
Boulder  Nonlinear  Systems  (BNS)  LC  SLM,  and  Santa  Barbara  Instruments  Group 
(SBIG)  camera.  The  reflective  SLM  has  an  array  size  of  7.68  mm  by  7.68  mm  con¬ 
sisting  of  512  by  512  active  pixels  and  a  2tt  phase  stroke.  The  reflected  wavefront 
distortion  is  A/ 15,  and  the  switching  frequency  is  30  Hz.  The  SBIG  charged  couple 
device  (CCD)  camera  consists  of  510  by  765  pixels  which  are  each  9  /im  by  9  //in. 

3.3  SLM  Characterization 

Per  design,  a  pixel  value  of  zero  ideally  results  in  maximum  birefringence  in  the 
device  while  the  strongest  field  and  minimum  birefringence  occur  with  a  pixel  value  of 
255.  The  intermediate  pixels  would  vary  linearly  from  the  minimum  to  the  maximum 
value;  however,  LC  molecular  rotation  does  not  vary  linearly  with  applied  voltage. 
Each  SLM  is  shipped  with  a  custom  lookup  table  (LUT)  to  give  an  approximate 
mapping  from  pixel  values  to  the  voltages  which  will  most  closely  result  in  a  linear 
phase  response.  The  phase  response  of  the  SLM  should  be  measured  to  ensure  the 
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appropriate  set  of  pixel  values  is  used  to  achieve  the  desired  optical  effect.  BNS 
guarantees  a  linear  phase  response  for  a  set  of  50-100  pixels. 

Fringe  patterns  are  created  using  the  Twyman-Green  interferometer  illustrated 
in  Figure  9.  The  HeNe  beam  passes  through  a  50/50  non-polarizing  beam  splitter 
(NPBS),  dividing  the  beam  into  one  that  illuminates  the  SLM  and  one  that  illumi¬ 
nates  a  reference  mirror.  The  SLM  and  mirror  are  aligned  so  that  they  are  nearly 
coplanar,  and  the  beams  interfere  to  form  fringe  patterns  on  the  charge-coupled  device 
(CCD)  camera. 


SLM 


SF 


HeNe 


-®  o 


NPBS 


I 

D 


Mirror 


CCD 


Figure  9.  The  Twyman-Green  interferometer  setup  uses  a  half-wave  plate,  beam  ex¬ 
panding  spatial  filter  (SF),  diaphragm  (D),  and  non-polarizing  beam  splitter  (NPBS). 
HeNe  light  from  the  SLM  and  mirror  legs  interfere  to  create  patterns  on  the  CCD. 


Two  techniques  for  measuring  the  phase  response  of  the  SLM  are  used  to  pin¬ 
point  the  linear  range  of  pixel  addressability.  In  the  first  case,  the  phase  response 
is  measured  by  loading  a  matrix  with  all  entries  consisting  of  one  constant  grayscale 
value  onto  a  large  portion  of  the  SLM  for  pixel  values  0-255  while  the  remainder  of 
the  SLM  is  left  to  act  as  a  mirror  as  depicted  in  Figure  10(a).  The  displacement  of 
each  fringe  pattern  for  the  modulated  portion  of  the  SLM  is  analyzed  with  respect 
to  the  unmodulated  portion.  In  the  second  test,  grayscale  images  consisting  of  al¬ 
ternating  columns  of  zero  value  and  values  ranging  from  0-255  are  loaded  onto  the 
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SLM  as  shown  in  Figure  10(c)  resulting  in  fringe  patterns  resembling  the  one  shown 
in  Figure  10(d). 


(a)  (b)  (c)  (d) 


Figure  10.  Calibration  fringes  (b)  and  (d)  result  from  (a)  constant  gray  scale  values  on 
a  portion  of  the  SLM,  and  (c)  gray  scale  images  with  constant  columns. 

Averaging  the  recorded  CCD  intensity  values  over  large  portions  of  the  resulting 
fringe  patterns  yields  sinusoidal  curves  from  which  phase  shift  can  be  calculated 
straightforwardly  with 


'=*{$)  <6> 
where  A  is  the  displacement  of  the  interference  fringe  peaks,  and  A  is  the  period  of 
the  interference  fringes  [13].  Figure  11  displays  the  phase  shift  at  each  of  the  256 
pixel  values  for  both  techniques  described,  normalized  so  that  the  zero  pixel  value 
corresponds  to  a  zero  phase  shift.  While  the  absolute  values  for  phase  shift  start  to 
diverge  as  pixel  values  get  higher  for  the  two  techniques,  they  have  a  similar  overall 
shapes. 

The  two  data  sets  are  averaged  and  smoothed  using  a  five-point  averaging  method 
to  capture  the  underlying  trend,  and  pixel  values  which  mostly  correspond  to  even 
increments  of  the  maximum  phase  delay  are  recorded  for  use  in  the  experiment.  As  a 
consequence,  more  of  the  pixel  values  used  come  from  steeper  segments  of  the  curve. 
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Figure  11.  Phase  response  data.  The  data  sets  corresponding  to  alternating  gray  scale 
values  and  the  application  of  a  constant  gray  scale  value  to  only  part  of  the  SLM  yield 
similarly  shaped  trends.  The  smoothed  average  is  used  to  pick  out  the  pixel  values 
giving  even  increments  of  the  available  phase  delay  range,  which  is  about  one  wave. 


3.4  Optical  Setup 

The  setup  used  by  Vellekoop  for  transmissive  inverse  diffusion,  illustrated  in  Fig¬ 
ure  12  involves  a  632.8  nm  HeNe  laser  beam  expanded  and  reflected  off  of  an  LC 
SLM.  Polarizers  and  quarter- wave  plates  are  used  for  operation  in  a  phase  mostly 
modulation  mode. 

The  primary  experimental  setup  used  for  reflective  inverse  diffusion  to  achieve  a 
focused  reflected  spot  was  designed  to  incorporate  principles  outlined  by  Vellekoop 
while  keeping  the  goal  of  aiding  in  indirect  photography  in  mind.  Like  the  experiments 
described  for  indirect  photography  and  transmissive  inverse  diffusion,  a  HeNe  laser  is 
used  as  the  light  source.  As  seen  in  Figure  13,  the  geometry  of  the  setup  is  similar 
to  that  of  transmissive  inverse  diffusion  with  the  exceptions  that  (1)  the  SLM  used 
is  phase  only  so  that  polarizers  are  not  needed,  (2)  the  microscope  objective  which 
collects  the  scatter  from  the  sample  is  removed  for  its  impracticality  with  respect  to 
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Figure  12.  Transmissive  inverse  diffusion  setup.  A  HeNe  beam  is  expanded  and  re¬ 
flected  off  an  LC  SLM.  Polarization  optics  select  a  phase-mostly  modulation  mode. 
The  SLM  is  imaged  onto  the  entrance  pupil  of  the  first  objective  and  focused  onto  the 
sample.  The  CCD  camera  images  the  transmitted  intensity  pattern. 

indirect  photography,  and  (3)  the  63x  objective  was  traded  out  for  a  lens  (L2)  with 
a  focal  length  of  1  meter.  Having  an  optical  component  so  close  to  either  side  of 
the  reflective  sample  would  render  reflective  inverse  diffusion  trivial  in  the  current 
context. 

The  30- mW  HeNe  beam  passes  through  a  half-wave  plate  to  orient  its  linear  po¬ 
larization  with  the  slow  axis  of  the  phase  only  SLM.  It  is  spatially  filtered,  expanded, 
collimated,  and  restricted  to  barely  underfill  the  SLM  before  reflecting  off  of  the 
mirror.  While  it  is  desirable  to  have  a  maximum  number  of  phase  modulating  pixels 
available  with  which  to  work,  overfilling  the  SLM  would  give  rise  to  uncontrolled  light 
reflecting  off  of  the  diffuse  surface.  The  beam  is  split  with  an  NPBS  which  dumps 
half  the  beam  and  transmits  the  other  half  for  reflection  off  of  the  cover  glass  of  the 
SLM.  The  beam  is  focused  onto  the  reflective  sample  whose  normal  is  oriented  45 
degrees  from  the  incoming  beam.  The  camera  is  placed  40  ±  0.5  cm  from  the  spot  on 
the  sample  at  45  degrees  on  the  other  side  of  the  surface  normal  so  that  it  is  centered 
on  the  specular  region  of  the  scattered  reflection.  All  the  reflective  surfaces,  detailed 
in  Section  3.8,  are  highly  scattering  samples. 
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Figure  13.  Reflective  inverse  diffusion  setup.  A  30  mW  HeNe  laser  beam  passes 
through  a  half-wave  plate,  is  expanded  (20x),  spatially  filtered  (SF),  and  collimated 
(LI),  and  is  restricted  to  the  size  of  the  SLM  with  a  diaphragm  (D).  It  is  reflected  off 
a  mirror  and  divided  with  a  non-polarizing  beam  splitter  (NPBS).  Half  reflects  off  the 
SLM,  split  again,  and  focused  (L2)  onto  the  sample.  The  sample,  which  is  oriented  at 
45  degrees  with  respect  to  both  the  incoming  beam  and  the  lensless  camera  (CCD), 
scatters  the  beam  onto  the  detector. 


3.5  Automated  Experimental  Control 

A  MATLAB  routine  was  developed  as  the  master  program  to  control  all  compo¬ 
nents  of  the  experimental  setup.  BNS  supplied  a  MATLAB  software  development 
kit  for  manipulating  the  SLM  in  the  form  of  a  graphical  user  interface  with  options 
for  preloading  a  set  of  phase  screens.  The  provided  code  was  modified  to  implement 
continuous  downloading  of  phase  screens  as  optimization  was  achieved.  The  SBIG 
CCD  was  controlled  in  MATLAB  through  calling  a  C-H-  executable  file  which  was 
written  and  compiled  in  Visual  Studio  2010.  In  a  single  iteration,  a  matrix  of  pixel 
values  was  loaded  onto  the  SLM,  the  camera  was  called  to  capture  an  image  with  a 
100-ms  exposure  time,  the  intensity  values  from  the  camera  were  read  into  MATLAB, 
the  values  from  the  target  focus  spot  were  summed,  and  the  resulting  target  intensity 
provided  an  indication  as  to  how  to  change  the  SLM  input  matrix  for  the  subsequent 
iteration.  Refer  to  Section  3.7  for  a  schematic  outlining  these  operations  and  the  time 
elapsed  for  each  operation. 
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3.6  Phase  Manipulation  Algorithms 


In  the  transmissive  inverse  diffusion  case,  there  is  a  unique  incident  wavefront 
which  makes  the  scattering  material  optimally  focus  light  to  a  specific  point.  The 
wavefront  is  phase  adjusted  on  the  scale  of  the  wavelength  of  light  and  cannot  be 
constructed  from  a  small  set  of  smooth  base  functions,  making  efficient  adaptive 
optics  algorithms  ineffective  [6].  As  a  result,  Vellekoop  implements  three  algorithms 
which  would  solve  for  the  phase  solution  by  computational  brute  force. 

The  first  is  a  stepwise  sequential  algorithm  which  varies  one  set  of  grouped  pixels 
at  a  time  to  find  and  save  the  pixel  value  which  leads  to  maximized  target  intensity  on 
the  detector,  then  loads  all  saved  pixel  values  onto  the  SLM  and  cycles  back  through 
for  another  iteration.  The  second,  depicted  in  Figure  14(c)  and  called  the  continuous 
sequential  algorithm,  is  the  same  as  the  first  except  that  the  saved  optimal  pixel  is 
loaded  onto  the  SLM  as  it  progresses  to  the  next  one.  The  final  is  a  partitioning 
algorithm  that  randomly  chooses  half  the  grouped  pixels,  varies  them  as  a  set,  and 
saves  the  optimal  pixel  value  in  the  entire  set;  it  is  illustrated  in  Figure  14(b)  [6]. 
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Figure  14.  Pixel  scheme,  (a)  SLM  pixels  which  touch  a  diameter  of  7.68  mm  when  the 
SLM  face  is  divided  into  a  grid  of  NxN  grouped  pixels  are  included  in  the  optimization 
scheme,  (b)  Half  the  pixels  are  randomly  chosen  and  optimized  as  a  set.  (c)  One  pixel 
is  varied  at  a  time  and  optimized  pixel  values  are  set  as  the  algorithm  runs. 


The  two  algorithms  depicted  in  Figure  14  are  used  in  reflective  inverse  diffusion 
as  a  result  of  their  reported  performance  results  in  transmissive  inverse  diffusion.  In 
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every  test  presented,  the  partitioning  algorithm  starts  achieving  enhancement  with 
the  fewest  iterations,  and  the  continuous  sequential  algorithm  achieves  the  highest 
overall  enhancement  by  the  end  of  each  run.  Figure  14(a)  illustrates  how  only  the 
SLM  pixels  for  a  given  number  of  addressable  segments,  which  the  diameter  of  the 
laser  touches,  are  varied  in  the  algorithm  while  all  remaining  pixels  in  the  corners 
are  left  at  a  zero  setting.  This  ensures  that  run  time  is  not  wasted  in  attempting  to 
optimize  segments  of  the  SLM  which  do  not  actually  affect  the  beam. 

The  pixel  radius  of  the  target  intensity  region  on  the  camera  is  set  to  the  approx¬ 
imate  size  of  a  speckle  for  every  test.  For  a  given  sample,  speckle  size  is  a  function 
of  the  number  of  modulated  segments  used.  A  smaller  segment  corresponds  to  a 
larger  number  of  total  segments  and  a  smaller  speckle  size  on  the  detector.  Chapter 
4  includes  an  illustration  of  how  speckle  size  varies  when  a  random  phase  map  of  size 
N=2  to  N=64  is  loaded  on  the  SLM  for  a  given  sample  and  optical  configuration. 
Especially  in  a  transition  from  one  set  of  segments  to  the  next,  the  speckle  size  will 
vary  as  optimization  is  taking  place.  To  account  for  variations  in  the  speckle  size 
for  a  given  number  of  segments,  more  weight  is  assigned  to  a  pixel  intensity  value 
located  at  the  center.  The  result  is  a  bulls  eye  pattern  of  weighted  target  intensities 
as  depicted  in  Figure  15. 


Figure  15.  Target  intensity  weighting.  The  target  area  on  the  detector  is  assigned  a 
weight  based  on  its  proximity  to  the  center.  The  most  highly  weighted  values  are  at 
the  center  (white)  and  decrease  toward  the  edge  (red). 
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Vellekoop  suggests  that  minimizing  the  effect  of  noise  can  be  accomplished  through 
coarse  pre-optimization  [14].  To  that  end,  the  algorithm  runs  with  small  sets  of  large 
grouped  pixels  before  decreasing  grouped  pixel  size  according  to  Figure  16.  Chapter 
4  opens  with  the  results  of  a  test  which  compare  the  enhancement  achieved  when  a 
phase  screen  of  a  particular  dimension  N  is  applied  with  and  without  optimizing  with 
the  preceding  smaller  dimension  N. 


Figure  16.  Grouped  pixel  dimensions.  For  a  given  reflector,  the  algorithm  runs  through 
increasing  grouped  pixel  dimension  size  TV,  from  N  =  2  to  N  =  32. 

3.7  System  Limitations 

The  most  notable  limitations  of  the  system  developed  to  demonstrate  reflective 
inverse  diffusion  center  on  the  speed  at  which  the  solution  space  is  searched  for  a 
phase  map  to  bring  the  scatter  to  a  focused  spot.  Figure  17  outlines  the  process 
by  which  the  target  area  is  optimized  and  the  time  it  takes  for  each  operation  in 
the  process.  Since  the  camera  is  opened  and  closed  to  take  a  single  frame  of  data 
between  each  iteration  of  the  SLM  being  loaded  with  a  particular  phase  screen,  the 
maximum  designed  refresh  rate  of  the  SLM  far  exceeds  the  refresh  rate  required  by 
the  algorithm  with  each  load. 

A  single  test  running  from  N  =  2  to  N  =  32  involves  the  optimization  of  a  total 
of  1160  individual  pixels.  Each  pixel  requires  an  image  be  taken  for  a  given  phase 
value.  If  the  phase  step  A(J)  is  set  to  7r/10,  each  pixel  must  run  through  21  values  for 
optimization.  The  test  thus  takes  a  little  over  ten  hours  to  run  through  a  set  for  a 
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Figure  17.  Algorithm  schematic  for  one  grouped  SLM  pixel.  To  take  a  single  measure¬ 
ment,  the  listed  algorithm  takes  approximately  1.5  seconds.  For  a  stepsize  A<1>  —  .  a 

grouped  pixel  is  optimized  in  around  32  seconds. 


single  sample.  Section  5.4  offers  suggestions  for  how  to  improve  upon  the  algorithm 
laid  out  here  in  order  to  achieve  the  desired  enhancement  with  significantly  less  run 
time. 

3.8  Reflection  Materials 

The  six  diffuse  reflectors  used  to  demonstrate  reflective  inverse  diffusion  include 
Spectralon,  a  dull  aluminum  sample,  a  slightly  polished  aluminum  sample,  Infragold, 
a  glass  cover  slip  covered  with  a  thin  layer  of  white  paint,  and  a  block  of  graphite. 
Figure  18  displays  photos  of  the  samples. 

Profilometer  measurements  are  taken  with  a  Tencor  Alpha  Step  1000  Surface 
Profiler  for  each  of  the  samples  used  in  reflective  inverse  diffusion.  A  single  scan 
lasts  8.5  minutes  while  covering  a  distance  on  the  surface  of  1  mm.  These  10  nm 
resolution  scans  are  included  in  Appendix  C,  and  their  corresponding  RMS  optical 
surface  roughnesses  are  calculated  with 

((z(x)-z)2dx  (7) 

where  L  is  the  length  of  the  scan  and  z(x)  is  the  surface  height  as  a  function  of 
position  x  [15]. 


25 


Figure  18.  Reflective  surfaces.  Reflective  inverse  diffusion  is  performed  on  and  com¬ 
pared  for  (a)  Spectralon,  (b)  unpolished  aluminum,  (c)  slightly  polished  aluminum,  (d) 
Infragold,  (e)  white  paint  on  a  glass  cover  slip,  and  (f)  graphite. 


The  target  focus  for  each  sample  is  aligned  with  the  center  of  the  specular  region 
of  reflection.  The  method  used  for  alignment  involves  starting  with  a  mirror  turned  to 
45  degrees  from  the  incident  beam  and  setting  up  a  fast  lens  to  focus  in  the  direction 
of  the  camera.  The  camera  is  translated  until  the  spot  is  focused  at  its  center.  The 
diffusely  reflecting  surface  is  inserted  and  rotated  to  45  degrees,  then  translated  until 
the  scattered  focus  is  centered  on  the  camera,  and  the  lens  is  removed.  The  process 
is  repeated  for  each  sample  because  the  mounts  and  reflector  thicknesses  vary  from 
sample  to  sample.  The  entire  setup  is  covered  with  black  cardboard  and  black  felt  to 
shield  all  optical  components  from  air  currents  and  external  reflective  surfaces. 

3.9  Nonmechanical  Beam  Steering 

Nonmechanical  beam  steering  in  the  context  of  reflective  inverse  diffusion  is  demon¬ 
strated  with  three  basic  experiments.  The  first  involves  consecutively  optimizing  to 
a  spot  in  each  corner  of  the  CCD  without  moving  anything  in  the  setup  using  the 
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white  paint  reflector.  The  second  and  third  require  the  translation  of  the  camera  from 
regions  of  specular  reflection  to  regions  of  diffuse  reflection.  Figure  19  illustrates  the 
simulation  of  the  scanning  of  a  playing  card  off  of  the  polished  aluminum  sample.  The 
specular  region  is  considered  the  starting  point  on  the  upper  left  corner  of  the  card. 
The  camera  is  translated  in  two  dimensions  so  that  the  distance  from  the  sample  to 
each  side  of  the  simulated  playing  card  is  fixed  at  40±0.25  cm. 


Figure  19.  Beam  steering  across  playing  card,  (a)  The  top  left  corner  is  set  to  the 
specular  reflection  region  of  the  scatter,  (b)  The  camera  is  translated  so  that  the 
detector  corresponds  to  the  corner  points  of  the  playing  card  while  keeping  the  distance 
between  the  reflector  and  the  CCD  at  40±0.25  cm. 


To  test  the  extent  of  the  ability  to  focus  the  diffuse  scatter  from  each  of  the 
samples,  the  camera  is  placed  40±0.25  cm  from  the  reflector  in  the  plane  of  incidence. 
This  corresponds  to  an  incident  angle  of  45±0.5  degrees  as  before,  and  a  scattering 
angle  of  0±0.5  degrees. 
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3.10  Speckle  Pattern  Decorrelation  Test 


Quantification  of  the  decorrelation  of  the  speckle  pattern  is  achieved  with  an  ex¬ 
periment  which  mimics  the  operations  involved  in  inverse  diffusion.  As  illustrated 
in  Figure  20,  the  algorithm  generates  and  loads  a  random  phase  screen  of  a  certain 
dimension  N.  The  phase  screen  is  saved  as  a  reference  phase  screen,  and  the  corre¬ 
sponding  image  of  the  speckle  pattern  on  the  detector  is  saved  as  the  reference  image. 
Random  phase  screens  are  generated  and  applied  to  the  SLM,  but  their  corresponding 
detector  images  are  not  saved.  The  algorithm  periodically  returns  to  the  reference 
phase  screen,  and  the  corresponding  image  is  saved.  The  time  at  which  the  reference 
phase  screen  is  reloaded  is  stored. 


Reference 
Phase  Screen 


Randomly  Generated 
Phase  Screens 


Repeated 
Phase  Screen 


Figure  20.  Speckle  decorrelation  test.  A  phase  map  of  a  certain  dimension  N  is  ran¬ 
domly  generated  and  saved  as  the  reference  SLM  while  its  corresponding  speckle  pat¬ 
tern  on  the  AD  is  saved  as  a  reference  image.  The  algorithm  continues  to  load  randomly 
generated  phase  maps  and  periodically  returns  to  the  reference  phase  map.  The  speckle 
patterns  which  correspond  to  the  reference  phase  map  are  recorded  for  comparison  with 
the  initial  recorded  speckle  pattern. 


The  correlation  coefficient  is  the  metric  used  for  assessing  the  similarity  between 
the  reference  image  and  subsequent  images  resulting  from  the  same  applied  phase 
screen.  Comparing  two  images  pixel  by  pixel  and  calculating  the  covariance  divided 
by  the  product  of  their  respective  standard  deviations,  the  coefficient  is  given  by 
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r  = 


(8) 


E„,  E„  (An.  ~  j)  (g.»n  ~  B) 
y (E,„ E„  (4™  -  3)2)  (e„ E„  (sm„  -  B)2) 

where  A  corresponds  to  the  reference  image  matrix  of  intensity  values,  B  corresponds 
to  the  subsequent  image  matrix,  the  subscript  m  denotes  the  pixel  row,  and  the 
subscript  n  denotes  the  pixel  column. 

3.11  Summary 

The  experiments  presented  in  this  thesis  aim  to  provide  a  proof  of  concept  for 
adapting  transmissive  inverse  diffusion  to  the  case  of  reflection,  particularly  in  geo¬ 
metric  setups  which  are  practical  within  the  context  of  indirect  photography.  The 
SLM  is  characterized  for  use  in  experiments  involving  optimization  to  achieve  en¬ 
hancement  in  the  specular  and  diffuse  regions  of  scatter  reflected  from  six  samples. 
The  same  algorithm  is  applied  on  transmission  for  comparison  to  the  literature  pre¬ 
sented  in  Chapter  2.  The  results  attained  using  the  methodology  given  in  this  chapter 
are  presented  and  analyzed  in  Chapter  4. 
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IV.  Results  and  Analysis 


4.1  Introduction 

The  results  outlined  in  this  chapter  include  the  reproduction  of  the  technique 
used  by  Vellekoop  in  focusing  light  transmitted  through  a  thin  white  paint  sample. 
These  measurements  serve  as  a  baseline  with  which  to  compare  results  achieved  on 
reflection  for  the  same  thin  paint  sample  along  with  the  brushed  aluminum,  sand¬ 
blasted  aluminum,  graphite,  spectralon,  and  infragold  reflectors.  For  every  reflector, 
enhancement  achieved  for  the  specular  region  is  further  compared  to  increasingly  dif¬ 
fuse  regions  of  scatter  in  order  to  evaluate  the  potential  for  steering  the  focused  spot. 
This  chapter  also  provides  insight  into  aspects  of  the  setup  and  technique  which  most 
affect  the  results  for  each  experiment. 

4.2  Optimization  Considerations 

Tests  are  performed  using  both  the  pre-optimization  scheme  outlined  in  Section  3.6 
and  without  pre-optimization  to  ensure  that  the  extra  time  used  to  optimize  smaller 
values  for  N  increases  rather  than  decreases  algorithm  efficiency.  Figure  21  displays 
the  results  for  a  set  of  runs  in  which  two  factors  are  varied  for  the  graphite  sample. 
In  the  test  denoted  by  the  solid  magenta  line,  the  phase  modulation  step  size  is 
the  smallest  to  achieve  greater  SLM  phase  resolution,  and  the  optimization  for  each 
dimension  occurs  after  being  pre-optimized  by  the  dimension  preceding  it.  In  the 
test  denoted  by  the  blue  dashed  line,  the  phase  modulation  step  size  is  reduced,  but 
it  still  undergoes  pre-optimization  for  each  dimension.  In  the  final  test  given  by  the 
black  dotted  line,  the  phase  modulation  step  size  is  reduced  by  the  same  amount  as 
in  the  preceding  test,  but  the  optimization  for  each  dimension  stands  alone  without 
pre-optimization. 
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Figure  21.  Pre-optimizing  the  SLM.  Intensity  profiles  cut  through  the  center  row  of 
the  AD  are  shown  for  the  optimization  of  reflection  off  of  graphite  for  pre-optimized 
runs  and  runs  that  are  not  pre-optimized.  Pre-optimizing  the  SLM  and  decreasing  the 
phase  delay  step  size  A<fi  show  higher  relative  intensity  of  the  optimized  spot.  The  plots 
come  from  the  final  intensity  profiles  for  the  (a)  N  =  4,  (b)  N  =  8,  (c)  N  =  16,  and  (d) 
N  =  18  optimization  schemes. 


Performance  of  the  algorithm  is  similar  for  all  cases  for  dimensions  TV  =  4  and 
N  =  8.  The  peak  intensities  start  diverging  for  N  =  16,  and  peak  intensity  actually 
decreases  from  the  N  =  16  to  the  N  =  32  stage  for  A<F  =  2tt / 5  without  pre¬ 
optimization.  As  a  consequence  of  these  results,  all  experiments  run  from  N  =  2  to 
N  —  32  consecutively  using  the  higher  resolution  phase  step.  An  added  advantage  of 
running  through  each  size  N  is  the  ability  to  pinpoint  the  phase  screen  dimension  for 
which  the  optimization  breaks  down  for  a  particular  sample  with  a  single  run. 

For  transmissive  inverse  diffusion,  Vellekoop  reported  that  the  target  intensity 
responds  sinusoidally  when  changing  the  phase  of  one  or  more  segments  on  the  SLM. 
The  wave  is  sampled  by  taking  ten  measurements  and  fitting  a  curve  which  exploits 
the  entire  range  of  available  phase  delay  [3].  For  reflective  inverse  diffusion,  the  current 
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setup  does  not  consistently  yield  a  sinusoidal  response  in  target  intensity  as  a  function 
of  phase  delay.  As  shown  in  Figure  22,  the  target  intensity  variation  from  infragold 
scatter  appears  roughly  sinusoidal  for  N  =  4  and  N  —  8,  but  varies  randomly  for  the 
other  stages.  All  samples  exhibit  similar  target  intensity  behavior. 


Multiple  of  tt/10  Phase  Delay  (rads) 


Figure  22.  Target  intensity  variation  with  phase  delay.  The  summed  target  intensity 
values  for  infragold  show  a  roughly  sinusoidal  pattern  for  N  =  4  and  N  =  8  while  the 
N  =  2,  N  =  16,  and  N  =  32  have  dissimilar  amplitudes  and  periods. 


As  a  result,  the  optimization  scheme  used  in  this  thesis  is  not  able  to  achieve  the 
resolution  used  by  Vellekoop  in  selecting  the  phase  delay  which  optimizes  a  particu¬ 
lar  pixel  on  the  SLM.  While  the  algorithm  takes  twice  as  many  measurements,  the 
the  algorithm  used  by  Vellekoop  further  fits  the  intensity  response  values  from  the 
detector  to  a  curve  given  by 

/  ($)  =  Ia  +  Ib  +  2^/TaIbcos  ($  —  $o)  (9) 

where  Ib  is  the  intensity  at  the  target  from  the  modulated  segments,  I  a  is  the  target 
intensity  caused  by  light  from  the  other  segments,  $  is  the  varied  phase,  and  $0  is 
the  unknown  optimal  phase  value  for  the  current  segment.  Vellekoop  is  therefore  able 
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to  assign  pixel  values  corresponding  to  the  entire  range  of  available  hardware-limited 
SLM  phase  delay. 

4.3  Speckle  Size 

Without  using  a  lens  to  collect  the  light  scattered  from  the  diffuse  surface,  control 
of  the  speckle  size  on  the  detector  is  limited  for  a  given  geometry.  Figure  23  illustrates 
how  with  varying  only  the  value  N  on  the  SLM,  the  speckle  size  changes  significantly. 
In  this  case,  the  sandblasted  aluminum  reflector  produced  the  interference  patterns 
shown  on  the  detectors  for  randomly  generated  phase  maps  with  N  =  2,  N  =  4, 
N  =  8,  N  =  16,  and  N  =  32.  All  samples  show  similar  variation  in  speckle  pattern 
size. 


Figure  23.  Speckle  size  varies  with  N.  Without  altering  the  geometry  of  the  setup, 
applying  random  phase  maps  with  (a)  N  =  2,  (b)  N  =  4,  (c)  N  =  8,  (d)  N  =  16,  and  (e) 
N  =  32  leads  to  significant  differences  in  the  speckle  size  of  the  interference  pattern  on 
the  detector.  The  patterns  displayed  here  are  from  the  dull  aluminum  reflector. 


As  mentioned  in  Section  2.4,  the  final  enhanced  spot  is  given  by  the  constructive 
interference  of  the  speckle  patterns  of  the  individually  phased  beams.  Section  refsec- 
SSA  relates  the  sizes  of  the  enhanced  spots  described  in  the  following  sections  to  the 
dimensions  of  the  phase  screens  used  to  formulate  them. 

4.4  Comparison  of  Reflective  Surfaces 

Table  1  outlines  the  performance  of  the  optimization  technique  for  each  of  the 
reflectors  along  with  its  measured  surface  roughness  Rrms  and  the  full  width  at  half 
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maximum  of  the  final  focused  spot.  The  spot  is  evaluated  using  Vellekoop’s  definition 
of  enhancement  in  Equation  (4).  It  is  calculated  by  picking  out  the  pixels  of  the 
detector  that  fall  within  the  full  width  at  half  maximum  (FWHM)  of  the  target  focus 
and  computing  the  collection’s  RMS  value  with  respect  to  the  mean  of  the  surrounding 
speckle  pattern.  The  enhancement  is  the  ratio  of  this  RMS  value  to  the  RMS  value  of 
the  surrounding  speckle  pattern  at  a  radius  of  three  times  the  FWHM  and  outward. 


Table  1.  Specular  Enhancement. 


Reflective  Samples 

Enhancement  (p#) 

Roughness  ( Rrms ) 

Final  FWHM 

Brushed  aluminum 

122.3 

1.5/im 

36±3pm 

Infragold 

89.9 

9.4pm 

38±3pm 

Sandblasted  aluminum 

67.7 

2.3pm 

38±3pm 

Graphite 

37.3 

3.5pm 

41±3pm 

White  paint 

36.8 

1.7pm 

41±3pm 

Spectralon 

13.8 

Unprofiled 

45±3pm 

Transmissive  Sample 

White  paint 

56.4 

1.7pm 

63±3pm 

The  enhancement  as  a  function  of  the  pixel  optimization  iteration  from  the  N  =  2 
through  the  N  =  32  cases  are  shown  in  Figure  24  with  the  N  =  16  and  N  =  32 
portions  highlighted  in  gray  and  yellow,  respectively.  Enhancement  for  every  reflector 
increases  most  rapidly  for  the  N=16  optimization  stage  and  either  decreases  or  levels 
off  during  the  beginning  of  the  N  =  32  stage.  Graphite  and  Spectralon  show  no 
improvement  during  the  N  =  32  stage,  and  white  paint  shows  little  improvement 
before  the  enhancement  decreases  again  to  a  value  slightly  above  its  enhancement  for 
N  =  16.  The  primary  factor  attributed  to  the  lack  of  enhancement  achieved  for  these 
three  samples  is  the  signal-to-noise  ratio  (SNR)  of  the  speckle  pattern  for  increasing 
SLM  dimension  N. 

A  second  factor  which  leads  to  further  decrease  in  enhancement  for  Spectralon 
as  the  algorithm  progresses  for  N  =  32  is  the  rate  at  which  the  speckle  pattern 
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Figure  24.  Sample  enhancement  comparison.  The  enhancement  r]R  as  defined  in  Sec¬ 
tion  2.4  is  plotted  as  a  function  of  the  optimized  pixel  iteration  for  SLM  dimension 
N  =  2  through  N  =  32.  The  TV  =  16  iterations  are  highlighted  in  gray  while  the  N  =  32 
iterations  are  highlighted  in  yellow. 


decorrelates.  Given  its  translucent  properties  mentioned  in  Chapter  2,  light  enters  the 
sample  and  undergoes  multiple  scattering  before  it  reemerges.  Recent  experiments 
involving  the  determination  of  the  penetration  depth  for  HeNe  light  in  samples  of 
Spectralon  have  produced  estimates  of  around  200/jin  [16].  This  value  is  around 
20  times  the  thickness  of  the  white  paint  sample,  indicating  that  the  shaped  HeNe 
beam  undergoes  many  more  scattering  events  upon  reflection  within  the  Spectralon 
sample  than  through  the  white  paint  sample  on  transmission.  A  specific  illustration 
of  the  decorrelation  of  the  Spectralon  speckle  pattern  on  reflection  and  the  white 
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paint  sample  speckle  pattern  on  transmission  as  a  function  of  time  is  shown  explicitly 
in  Section  4.7. 

The  brushed  aluminum  reflector  yielded  the  highest  enhancement  of  the  set.  Fig¬ 
ure  25(a)  shows  the  intensity  ratios  between  the  final  images  for  each  phase  map 
pixel  dimension  N  and  the  average  starting  speckle  pattern  for  that  N .  Random 
phase  maps  for  a  particular  size  N  are  applied  to  determine  a  starting  intensity  value 
to  which  the  images  and  plots  are  normalized.  The  speckle  patterns  are  displayed 
using  a  logarithmic  scale  in  order  to  capture  the  intensity  distribution  around  the 
focused  spot.  See  Appendix  D  for  similar  displays  of  the  other  reflecting  surfaces. 

Figure  25(b)  shows  the  same  information  for  the  optimization  process  for  the 
white  paint  sample  on  transmission.  The  speckle  size  remains  fairly  constant  as  the 
SLM  dimension  size  N  is  increased.  The  optimization  algorithm  fails  to  increase 
the  enhancement  for  the  sample  for  the  N  =  32  stage  due  to  the  decorrelation  of 
the  speckle  pattern.  The  N  =  32  curve  shows  a  smearing  out  of  the  optimized 
spot  in  spite  of  a  consistent  background  speckle  size.  This  can  be  thought  of  as 
the  individual  speckle  patterns  each  decorrelating  to  smear  out  the  effect  of  their 
constructive  interference  at  the  target. 

4.5  Spot  Size  Analysis 

The  varying  speckle  size  observed  in  Section  4.3  leads  to  corresponding  variations 
in  the  final  focused  spot  size  for  each  SLM  dimension  N.  The  relationships  which 
govern  these  variations  come  from  considering  the  light  propagating  from  a  given 
pixel  as  an  individual  beam  and  applying  Fraunhofer  diffraction  theory.  Figure  26 
illustrates  the  relevant  quantities. 

The  spot  diameters  q\  and  q-2  on  the  reflective  surface  and  on  the  detector,  respec¬ 
tively,  go  as 
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Figure  25.  Sample  relative  intensity  comparison.  The  top  half  of  each  block  displays 
the  interference  pattern  from  a  portion  of  the  detector  center  for  phase  screens  of 
N  =  4,  N  =  8,  N  =  16,  and  N  =  32.  The  pattern  is  normalized  to  the  average  starting 
intensities  and  plotted  on  a  logarithmic  scale.  An  image  of  the  corresponding  phase 
screen  applied  is  shown  with  a  phase  delay  of  zero  corresponding  to  black  and  a  phase 
delay  of  2i r  corresponding  to  white.  The  bottom  half  of  each  block  shows  the  center 
row  of  intensity  values  recorded  from  the  detector  on  the  final  image  of  each  SLM 
dimension  optimization.  This  data  is  displayed  for  (a)  brushed  aluminum  on  reflection 
and  (b)  white  paint  on  transmission. 
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CCD 


Figure  26.  Spot  size  relationships.  The  spot  size  on  the  sample  q\  is  a  function  of  the 
grouped  pixel  size  dpix.  This  in  turn  affects  the  size  of  the  enhanced  spot  q2  recorded 
on  the  detector. 
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Substituting  for  dPiX  and  simplifying,  the  final  relationship  becomes 
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Considering  that  R\  ~  /,  the  lens  focal  length,  geometry  of  the  setup,  and  the  pixel 
dimension  are  the  primary  factors  in  predicting  the  spot  size  of  the  enhanced  spot  on 
the  detector. 

In  order  to  extract  the  boundaries  of  the  focused  spot  from  the  background  of 
random  speckle  and  noise,  a  1-dimensional  autocorrelation  is  applied  to  the  center 
row  of  intensity  values  from  the  AD.  The  original  intensity  distribution  as  a  function 
of  x  is  displaced  by  an  amount  e.  Displacing  the  function  a  pixel  at  a  time,  the  area 
under  the  product  of  the  displaced  and  undisplaced  versions  of  the  distribution  give 
the  degree  of  overlap  in  accordance  with 


f  (x  +  e)  f*  ( x )  dx. 


(13) 


cll 


The  function  f*  (x)  denotes  the  complex  conjugate  of  f  (x). 

Figure  27  displays  the  autocorrelations  performed  for  every  reflective  sample  for 
the  final  optimization  iteration  corresponding  to  N  =  8,  N  =  16,  and  N  =  32. 
The  distance  between  e  =  0  and  the  point  at  which  the  autocorrelation  values  start 
leveling  off  into  a  baseline  noise  value  gives  the  full  width  of  the  enhanced  spot.  The 
vertical  lines  point  to  the  values  of  the  predicted  spot  size  (red  center  lines)  and 
the  uncertainty  bounds  (cyan  and  purple  outer  lines).  The  autocorrelations  show 
similar  spot  sizes  for  each  pixel  dimension  for  all  six  reflective  samples  and  are  in 
good  agreement  with  the  predicted  values  of  384  /un  for  N  =  8,  192  /an  for  N  =  16, 
and  96  /im  for  IV  =  32. 


Figure  27.  Signal  autocorrelation.  The  autocorrelation  is  shown  for  each  sample  for 
SLM  dimensions  N  =  8  through  N  =  32.  The  vertical  bands  denote  the  predicted  value 
(red  center  line)  with  uncertainty  bounds  (cyan  and  purple  lines)  of  the  spot  diameter. 
The  point  at  which  the  autocorrelation  levels  out  corresponds  to  the  edge  of  the  focused 
spot. 


4.6  Nonmechanical  Beam  Steering 

As  mentioned  in  Chapter  2,  nonmechanical  beam  steering  is  typically  thought  of 
as  a  technique  to  steer  a  coherent  beam  effectively  using  a  phase  screen  as  a  diffraction 
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grating.  In  the  context  of  reflective  inverse  diffusion,  nonmechanical  beam  steering  is 
equivalently  the  ability  to  choose  a  displaced  target  area  to  overlay  speckle  patterns 
of  individually  modulated  segments  of  the  beam  and  achieve  a  focused  spot.  The 
first  demonstration  of  beam  steering  involves  choosing  the  four  corners  of  the  CCD 
to  optimize  reflection  from  the  white  paint  sample  using  N  =  2  to  IV  =  16.  The 
optimization  algorithm  for  each  spot  is  the  same  as  that  outlined  in  Section  3.6, 
except  that  it  does  not  progress  through  the  N  =  32  stage,  and  the  algorithm  runs 
four  times  consecutively  with  target  areas  defined  in  one  of  the  four  corners  instead  of 
the  center.  Figure  28  shows  the  four  final  iteration  images  using  the  same  logarithmic 
scaling  and  normalization  to  the  average  starting  intensity  as  is  used  in  Figure  25. 


Figure  28.  Beam  steering  to  CCD  corners.  Continuous  runs  from  iV  =  2  to  iV  =  16  are 
performed  on  the  white  paint  sample  for  the  four  corners  of  the  CCD  to  illustrate  how 
the  detector  can  be  scanned  without  moving  anything  in  the  setup. 

As  expected,  all  spots  show  similar  peak  intensity  values.  They  are  separated 
lengthwise  by  6.9  mm,  4.6  mm  vertically,  and  8.3  mm  diagonally.  This  corresponds 
to  a  steer  angle  of  31.3  mrad  with  respect  to  the  scattering  surface.  The  peak  intensity 
of  the  least  enhanced  spot  on  the  CCD  is  within  13.5%  of  that  of  the  most  highly 
enhanced  spot. 
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Steering  across  the  face  of  a  playing  card  requires  a  beam  displacement  of  88 
mm  vertically,  63  mm  horizontally,  and  108.2  mm  diagonally.  Figure  29  displays  the 
intensity  profile  slices  corresponding  to  spots  at  each  of  these  corners  as  described  in 
Section  3.9  for  polished  aluminum.  The  plots  are  arranged  in  the  figure  according  to 
their  locations  on  a  the  face  of  a  playing  card;  i.e.,  the  top  left  plot  corresponds  to  the 
top  left  corner  of  the  card,  etc.  Enhancement  decreases  by  58.6%  with  the  diagonal 
distance  from  the  top  left  specular  region  of  the  card  to  the  bottom  right  region. 
The  reduction  in  enhancement  is  smaller  than  the  plotted  reduction  in  peak  intensity 
because  the  enhancement  is  the  ratio  of  the  RMS  focus  intensity  to  the  surrounding 
speckle  pattern,  and  the  RMS  intensity  of  the  surrounding  speckle  pattern  for  the 
specular  region  of  the  playing  card  is  greater  than  for  the  diffuse  region. 

To  compare  the  extent  of  beam  steering  capabilities  and  to  allow  for  comparison 
across  all  samples,  the  CCD  is  moved  to  fall  in  line  with  the  surface  normal,  cor¬ 
responding  to  an  incident  angle  of  45±.5  degrees  and  a  scattering  angle  of  0±0.5 
degrees.  Figure  30  shows  the  enhancement  as  a  function  of  pixel  optimization  itera¬ 
tion  for  each  of  the  six  samples.  Both  the  specular  reflection  enhancement  and  the 
diffuse  reflection  enhancement  are  plotted. 

For  all  samples,  the  test  runs  through  IV  =  2  to  N  =  16.  Polished  aluminum,  the 
reflector  which  yielded  the  highest  specular  reflection  enhancment,  achieved  negligible 
enhancement  in  the  more  diffuse  region  of  its  scatter  due  to  its  overwhelming  concen¬ 
tration  of  specular  reflection.  The  sandblasted  aluminum  and  infragold  enhancements 
in  the  direction  of  the  surface  normal  reached  45.40%  and  40.04%  of  their  specular 
enhancements,  respectively. 

Spectralon,  a  highly  Lambertian,  reflective  sample  with  a  short  persistence  time 
performed  nearly  equivalently  for  both  specular  and  diffuse  enhancement.  Graphite, 
a  good  absorber  in  the  visible  spectrum,  showed  negligible  enhancement  in  the  diffuse 
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Figure  29.  Beam  steering  along  dimensions  of  a  playing  card.  The  polished  aluminum 
reflector  is  used  to  illustrate  the  decrease  in  enhancement  for  spots  corresponding  to 
the  corners  of  a  playing  card.  The  top  left  corner,  the  optimization  target  which  lies 
in  the  specular  region  of  reflection,  achieves  the  highest  enhancement.  Enhancement 
falls  off  with  distance  from  the  specular  region. 


region.  Noted  earlier  for  its  loss  due  to  transmission,  the  white  paint  sample  showed 
little  enhancement  in  the  diffuse  region,  as  well. 

For  high  displacement  beam  steering,  the  samples  which  are  both  strongly  scat¬ 
tering  and  highly  reflective  tend  to  achieve  enhancements  more  comparable  to  their 
respective  specular  enhancements.  The  use  of  nonmechanical  beam  steering  is  thus 
limited  by  the  intensity  of  light  reflected  at  various  scatter  angles.  In  the  context  of 
indirect  photography,  for  which  an  object  is  potentially  scanned  with  a  focused  beam 
using  reflective  inverse  diffusion,  the  scattering  properties  of  diffuse  materials  become 
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Figure  30.  Enhancement  comparison  of  specular  to  diffuse  reflection.  Enhancement 
is  plotted  as  a  function  of  the  optimization  iteration  for  each  sample  with  the  camera 
placed  in  the  specular  region  of  reflection  (blue)  and  at  a  scatter  angle  of  0±0.5  degrees 
(green). 


an  aid  rather  than  a  hindrance  for  objects  extending  outside  of  the  specular  region  of 
scatter. 

4.7  Speckle  Pattern  Decorrelation 

Further  exploration  of  the  speckle  pattern  decorrelation  shows  trends  which  vary 
from  sample  to  sample.  Correlation  coefficients  between  a  reference  image  of  the 
speckle  pattern  for  a  particular  sample  to  subsequent  images  of  the  speckle  pattern 
are  computed  as  outlined  in  Section  3.10.  Plotting  these  coefficients  as  a  function  of 
the  time  which  passes  between  the  reference  image  and  the  comparison  images  being 
taken  gives  rise  to  Figure  31. 
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Figure  31.  Speckle  decorrelation.  Images  taken  for  a  particular  phase  screen  are  com¬ 
pared  to  the  initial  reference  image  using  the  correlation  coefficient,  and  the  coefficients 
are  plotted  as  a  function  of  time. 


It  should  be  noted  that  the  displayed  results  of  the  decorrelation  test  are  depen¬ 
dent  on  both  the  properties  of  the  material  and  the  environmental  conditions  for  each 
individual  test.  The  oscillatory  nature  of  some  of  the  data  suggests  that  certain  peri¬ 
odic  perturbations  exist  within  the  setup,  e.g.  fluctuations  in  detector  fan  speed,  laser 
output  power,  or  laboratory  HVAC.  Section  5.4.2  offers  suggestions  for  incorporating 
speckle  decorrelation  data  to  characterize  a  particular  experiment. 

4.8  Summary 

Starting  intensity  on  the  detector  and  the  sample  stability  are  the  main  factors 
determining  the  level  of  enhancement  that  can  be  achieved  through  this  optimization 
process.  This  implies  that  the  reflectance  properties  of  a  particular  sample  determine 
the  potential  for  enhancement  in  the  specular  direction  and  for  beam  steering  ap- 


44 


plications.  Sample  stability  given  the  current  laboratory  conditions  can  be  assessed 
using  the  correlation  coefficient  when  comparing  images  of  the  speckle  pattern  for 
a  given  SLM  phase  screen.  Geometric  considerations  such  as  speckle  size  become  a 
factor  when  lenses  with  longer  focal  lengths  are  used  to  focus  the  shaped  beam  onto 
the  sample.  For  real-world  application  of  this  technique  within  indirect  photography, 
the  ultimate  desire  is  to  further  distance  the  optical  system  from  the  reflecting  surface 
and  achieve  enhancement  in  a  variety  of  environmental  conditions. 
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V.  Discussion 


5.1  Introduction 

The  research  efforts  detailed  in  this  document  are  aimed  at  providing  a  proof 
of  concept  for  the  focusing  of  light  scattered  from  diffuse  reflectors  using  an  SLM 
to  assign  phase  delay  to  the  impinging  beam.  All  experiments  are  performed  in 
setups  which  are  practical  within  the  context  of  indirect  photography.  The  four  main 
objectives  which  were  introduced  in  Chapter  1  include  this  demonstration  of  a  proof  of 
concept,  the  nonmechanical  steering  of  the  focused  spot,  the  identification  of  certain 
reflective  properties  which  affect  the  results,  and  an  outline  detailing  the  progression 
required  to  advance  this  technique  for  implementation  in  indirect  photography.  This 
chapter  serves  to  summarize  the  findings  and  resulting  lessons  learned  and  sketch  out 
steps  to  be  taken  toward  incorporating  this  technique  into  indirect  photography. 

5.2  Results  Summary 

Using  the  algorithm  developed  in  Section  3.6  and  the  enhancement  definition 
established  in  Section  2.4,  the  enhancement  achieved  on  transmission  for  the  white 
paint  sample  was  56.4.  The  more  specular  samples  achieved  enhancements  between 
67.7  and  122.3  in  the  specular  direction.  The  more  diffuse  samples  achieved  lower 
enhancements  in  the  specular  direction,  but  higher  enhancements  in  the  direction  of 
the  surface  normal  relative  to  their  corresponding  specular  direction  enhancements 
for  a  beam  incident  at  45°. 

5.3  Conclusions  Drawn  from  Results 

The  demonstration  of  an  SLM  introducing  phase  delay  patterns  which  inversely 
diffuse  a  HeNe  beam  to  correct  for  diffuse  scattering  off  of  a  reflective  surface  was 


46 


accomplished.  High  specular  reflectance  and  sample  stability  contribute  most  signifi¬ 
cantly  to  achieving  enhancement  in  the  specular  region,  while  diffuse  scattering  angle 
reflectance  is  significant  in  achieving  comparable  diffuse  region  enhancement.  There 
is  still  a  great  deal  of  work  to  be  done  to  implement  this  technique  for  reducing  the 
problem  of  indirect  photography  to  an  application  of  dual  photography,  as  outlined 
in  Section  5.4. 

5.4  Recommendations  for  Future  Work 

The  recommendations  for  future  work  fall  into  three  broad  categories:  (1)  ways  to 
improve  the  current  setup  to  achieve  higher  enhancement  values,  (2)  additional  work 
required  to  more  fully  characterize  the  physics  underlying  the  technique,  and  (3)  the 
steps  required  to  incorporate  reflective  inverse  diffusion  into  indirect  photography. 

5.4.1  Algorithm  Improvement. 

Any  change  to  the  algorithm  which  searches  the  solution  space  and  pinpoints 
a  phase  solution  more  quickly  will  incite  greater  enhancement  in  the  focused  spot. 
The  time  it  takes  to  achieve  a  solution  will  always  be  at  odds  with  the  persistence 
time  of  the  reflected  speckle  pattern.  Calling  the  camera  to  take  a  single  image 
and  reading  target  intensity  values  into  MATLAB  from  any  of  the  SBIG  camera’s 
default  photo  formats  take  an  exorbitant  amount  of  time.  The  camera  does  not 
remain  open  while  new  SLM  phase  screens  are  generated,  thus  none  of  the  operations 
are  synchronized  for  efficiency.  Given  the  same  lab  equipment,  potential  methods  for 
better  integrating  the  controlling  software  include  but  are  not  limited  to:  (1)  building 
an  interface  using  a  MATLAB  adaptor  kit  which  is  supplied  with  the  MATLAB 
Image  Acquisition  Toolbox.  At  a  minimum,  this  requires  a  substantial  knowledge 
of  the  C++  programming  languages.  (2)  Create  a  MEX  interface  with  information 
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provided  at  the  Mathworks,  Inc.  website.  (3)  Work  with  the  software  development 
consultants  at  AFIT  to  write  the  interface. 

Establishing  a  routine  which  minimizes  the  time  to  focus  to  a  sufficiently  intense 
spot  for  scanning  the  hidden  object  is  desirable.  Other  options  for  decreasing  run  time 
include  implementing  adaptive  algorithm  techniques  for  approaching  an  improved 
solution.  As  an  example,  for  segments  which  do  not  contribute  much  to  improving 
the  target  intensity,  the  algorithm  could  break  up  regions  of  the  phase  map  into 
smaller  segments.  Figure  32  shows  a  16x16  phase  screen  with  a  region  broken  into 
segments  corresponding  to  a  32x32  phase  screen. 


Figure  32.  Adaptive  algorithm.  Given  certain  indicators  in  the  target  intensity  for  a  set 
of  grouped  pixels,  break  the  pixel  into  subsets  of  smaller  grouped  pixels  to  efficiently 
increase  enhancement. 

5.4.2  Technique  Improvement. 

The  calibration  technique  used  to  characterize  the  phase  modulation  of  the  SLM 
assumed  a  flat  phase  response  across  the  SLM.  In  an  effort  to  glean  a  set  of  pixel  values 
which  would  best  cover  the  range  of  phase  delay,  sets  of  fringe  patterns  corresponding 
to  opposing  gray  scales  were  compared.  To  ensure  that  each  and  every  pixel  most 
closely  covers  the  range  of  phase  delay  in  even  increments,  one  can  further  calibrate 
the  SLM  static  aberration.  Most  often  and  notably  caused  by  a  backplane  which 
is  curved  from  manufacturing  processes,  static  aberration  can  be  tested  for  using 
interferogram  analysis,  a  commercial  interferometer,  or  parameterized  phase  retrieval 
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[17,  18,  13].  Interferogram  analysis  can  be  done  with  a  Twyman-Green  interferometer 
setup  similar  to  the  one  used  for  fringe-shift  analysis.  This  analysis  would  result  in 
a  LUT  that  is  different  for  each  pixel.  To  reduce  the  size  of  the  LUT  associated 
with  the  SLM  for  pixel-by-pixel  conversion  factors,  the  impinging  beam  size  could  be 
reduced  to  cover  flatter  regions  of  the  SLM. 

The  performance  of  the  enhancement  process  is  highly  dependent  on  the  starting 
speckle  intensity.  The  current  setup  is  designed  to  resemble  the  setup  used  for  inverse 
diffusion  on  transmission  for  easy  comparison;  however,  half  of  the  beam  is  dumped 
on  the  first  pass  through  the  NPBS  onto  the  SLM,  and  more  is  lost  upon  its  second 
reflection  off  the  NPBS.  Implementing  an  off-axis  setup  will  ensure  that  the  more 
power  is  delivered  to  the  surface  reflector  for  a  given  laser.  Power  considerations  will 
become  even  more  significant  as  the  optical  system  is  pushed  back  farther  from  the 
reflecting  surface. 

The  decorrelation  data  presented  in  Section  4.7  was  taken  to  characterize  how  the 
speckle  pattern  for  a  particular  phase  map  decorrelates  for  a  given  sample  and  set  of 
environmental  conditions.  The  operations  of  an  enhancement  experiment  are  mim¬ 
icked,  suggesting  that  if  an  enhancement  experiment  had  been  performed  instead,  it 
would  have  suffered  from  the  amount  of  decorrelation  recorded.  In  order  to  disentan¬ 
gle  the  effect  of  the  speckle  decorrelation  due  to  sample  and  environmental  stability 
from  the  effect  of  changing  other  factors  of  interest,  one  could  incorporate  the  decor¬ 
relation  measurement  technique  directly  into  the  enhancement  technique.  Defining 
a  particular  random  phase  map  as  a  reference  and  periodically  loading  it  onto  the 
SLM  throughout  the  enhancement  process  would  provide  the  information  necessary 
to  establish  a  decorrelation  metric  for  each  test. 
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5.4.3  Reflective  Inverse  Diffusion  Characterization. 


Additional  work  must  be  done  to  pin  down  the  contribution  of  various  factors 
affecting  the  optimization  process.  All  reflective  inverse  diffusion  experiments  were 
performed  using  a  lens  with  a  1-meter  focal  length  to  focus  the  phase  modulated 
beam  onto  the  reflective  sample.  To  glean  a  relationship  between  the  spot  size  on  the 
sample  and  the  resulting  enhancement,  lenses  with  differing  focal  lengths  should  be 
used  and  their  results  compared. 

Another  important  relationship  may  be  the  effect  of  varying  starting  intensity  on 
the  detector  for  a  given  sample  within  a  given  geometry.  In  general,  samples  with 
higher  initial  speckle  pattern  intensities  achieved  better  results;  however,  exploring 
the  relationship  between  starting  intensity  and  resulting  enhancement  for  a  particular 
sample  or  among  various  samples  may  provide  insight  into  enhancment  limits  imposed 
by  certain  material  properties.  Furthermore,  a  connection  has  been  made  between  the 
BRDF  and  the  resulting  enhancement  achieved  for  the  samples  involved.  A  precise 
relationship  between  a  sample’s  BRDF  and  enhancement  can  be  made  by  comparing 
enhancement  for  a  range  of  angles  to  the  scatter  data  provided  by  a  Complete  Angle 
Scatter  Instrument  (CASI). 

5.4.4  Implementation  in  Indirect  Photography. 

As  mentioned  in  the  introduction,  there  is  a  progression  of  advancements  in  this 
research  that  must  be  made  before  reflective  inverse  diffusion  could  be  incorporated 
to  aid  in  indirect  photography.  The  aim  of  indirect  photography  is  to  image  objects 
hidden  from  both  the  camera  and  the  light  source;  therefore,  the  SLM  and  detector 
which  provides  feedback  to  the  SLM  should  not  have  access  to  a  direct  line  of  sight  to 
the  object.  After  improving  the  algorithm  to  achieve  a  higher  intensity  focused  spot 
on  the  camera,  which  currently  takes  the  place  of  the  object  to  be  imaged,  a  second 
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detector  must  be  introduced  to  monitor  interference  patterns  as  the  SLM  modulates 
the  beam.  Eventually,  the  detector  should  be  placed  in  close  proximity  with  the  other 
components  of  the  indirect  photography  system  while  giving  feedback  to  the  SLM  to 
focus  onto  spots  on  the  object.  An  intermediate  step  could  involve  dividing  the 
camera  into  two  regions  and  focusing  to  a  spot  on  one  side  while  providing  feedback 
to  the  SLM  from  the  other. 

5.5  Final  Comments 

The  research  detailed  in  this  thesis  completes  the  first  step  in  reducing  the  com¬ 
plexity  of  the  imaging  process  of  indirect  photography  to  the  process  used  in  dual 
photography.  It  establishes  a  proof  of  concept  for  focusing  the  light  scattered  off  of 
diffuse  surfaces  using  phase  modulation  and  steering  that  focus.  There  is  a  signifi¬ 
cant  amount  of  work  to  be  done  to  advance  this  technique  for  realistic  use  in  indirect 
photography,  including  improving  the  current  technique,  establishing  precise  relation¬ 
ships  between  enhancement  and  experimental  factors,  and  eventually  collocating  the 
detector  providing  feedback  to  the  SLM  with  the  rest  of  the  imaging  system. 
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Appendix  A.  Reflective  Inverse  Diffusion  MATLAB 

Function 


Appendix  A  provides  the  main  function  which  controls  the  SLM  and  the  camera. 
The  subroutines  are  listed  separately  in  Appendix  B  so  that  the  root  function  serves 
as  a  clear  outline  of  the  operations  performed  for  the  experiment. 


function  ExperimentNMat rix  (Folder , PixelFile, NMatrix, PixRadius ) 

%  Written  by  Jessica  Schafer,  2Lt,  AFIT/ENP 
%  Updated:  15  October,  2011 

%  All  functions  starting  with  ”BNS_"  are  taken  from  the  MATLAB  SDK  provided 
%  by  BNS  with  minor  alterations  for  the  current  context.  NMatrix  is  a 
%  vector  of  phase  map  dimensions.  The  most  common  input  for  this  thesis  is 
%  [2  4  8  16  32] .  PixRadius  is  determined  by  sample  speckle  size.  PixelFile 
%  is  the  set  of  pixel  values  which  will  give  the  desired  phase  delay. 

%%  Initialize  SLM. 

BNS  .Open  SLM; 

BNS_SetPower  (1)  ; 

ImageFrame  =0;  %  Set  PCI  card  memory  slot  value 

%%  Read  in  files  which  will  be  used  in  routine. 

LUTMap  =  BNS_ReadLUTFile ( 'C:\BNSMatlabSDK\LUT_Files\linear.LUT' ) ; 

PixValues  =  dlmread (PixelFile ) ; 

%%  Initialize  parameters. 

WinningMat rix=zeros (510,  765)  ; 

ImageMatrix  =  round (255*rand (512 )) ; 

Int=zeros (1, length (PixValues) ) ; 

Count  =  0; 

Compact ImageMatrix  =  round ( 255*rand (NMatrix ( 1 ))) ; 
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%%  Record  initial  captures  for  each  N  value 
RecordAvgStart Intensities  (NMatrix, LUTMap, Folder) 

%%  Vary  phase  pixel  by  pixel 

for  h  =  1 : length (NMatrix) 

N  =  NMatrix (h) ; 
n  =  num2str (N) ; 

H  =  num2str (h)  ; 
for  i=l:N 

I=num2str (i, '  % . 2d '  )  ; 
for  j=l : N 

%%  Only  vary  pixel  if  it  meets  condition  determining  if  it 

%%  falls  within  a  circle  of  radius  7.62  mm. 

if  (  (i<N/2&& j<N/2) && ( sqrt ( (N/2-i)  ~2+ (N/2-j) ~ 2 )  <N / 2  )  )  1 1  .  .  . 

( (i>N/2&& j<N/2) && (sqrt ( (N/2- (i-1) ) ~2+ (N/2-j) *2)  <N/ 2  )  )  1 1 . . . 

( (i<N/2&& j>N/2) && (sqrt ( (N/2-i) ~2+ (N/2- ( j-1) ) ~2) <N/  2  )  )  1 1  . . . 

( (i>N/2&& j>N/2) && (sqrt ( (N/2- (i-1) ) ~2+ (N/2- (j-1) ) ~2)<N/2) ) 

for  k=l : length (PixValues )  ; 

[Compact ImageMatrix, ImageMatrix] = . . . 

ChangePhaseMap ( ImageMatrix, Compact ImageMatrix, i, . . . 
j , PixValues (k) , N) ; 

BNS_LoadImageFrame (ImageFrame, ImageMatrixf  LUTMap) ; 
BNS_SendImageFrameToSLM (ImageFrame)  ; 

ImageFrame=mod ( ImageFrame  +  1 ,  26); 

Capturelmage ; 

CaptureMat rix=ReadSaveImage ( ' Capture . sbig ’ ) ; 

Int (k) =CalculateIntensityProf ile (CaptureMat rix, . . . 
PixRadius (h) ) ; 

%%  Save  the  pixel  value  and  corresponding  image  which 
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%%  gives  the  highest  intensity  according  to  the  target 

%%  function 

if  Int (k) ==max ( Int ) 

WinningMat rix  =  CaptureMat rix; 

end 


J=num2str ( j ,  1  % . 2d ’  )  ; 

if  k==length (PixValues ) 

Count  =  Count  +  In¬ 
compact  ImageMatrix  =  ... 

DeterminePixelValue (PixValues, Int, i, j, . . . 
Compact ImageMatrix) ; 
dlmwrite (strcat . . . 

(’C:\BNSMatlabSDK\Schafer  SLM  Cont rol\ ' , . . . 
Folder, ’ \ Opt Image '  ,  I ,  J,  n,  H, ’ . txt ' ) , . . . 
WinningMat rix,  ' delimiter '  ,  1  \t ' ) ; 

imagesc (WinningMat rix) ;  colormap ( ’ hot ' ) ; 
t it le  ( strcat ( ' Step  1 , I , J, n, H) ) ; 
dlmwrite (strcat . . . 

('C:\BNSMatlabSDK\Schafer  SLM  Cont rol\ ’ , . . . 
Folder, ' \ Tar get Intensity '  ,  I ,  J,  n,  H, . . . 

'  . txt ' ) , Int ,  ' delimiter '  ,  ' \t  ’  ) ; 

end 

end 

end 

end 


end 

%%  Save  the  phase  screen  which  optimizes  the  target  intensity, 
dlmwrite ( strcat ( 1 C : \BNSMat labSDK\ Schafer  SLM  Control\ ’, Folder , . . . 

1 \PhaseScreen ' , n, H,  ’  . txt ' ) , Compact ImageMatrix,  ' delimiter  1 ,  ' \t ' ) ; 

%%  Expand  compact  image  matrix  of  the  previous  size  N  to  the  next  size  N. 
if  hylength (NMatrix) 
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Compact ImChange  =  zeros (NMatrix (h+1 )) ; 
for  col  =  1 : N 

for  row  =  1 : N 

Compact ImChange ( (2*row—  1)  :  (2*row) ,  (2*col— 1)  :  (2*col) )  .  .  . 
=Compact ImageMatrix (row, col)  ; 

end 

end 

Compact ImageMatrix  =  Compact ImChange; 

end 

end 

%%  Display  last  image. 

imagesc (CaptureMat rix) ;  colormap ( ’ hot ’ ) ; 

figure;  imagesc (Compact ImageMatrix) ;  colormap (' gray  1 ) ; 

BN  S  _C  1  o  s  e  S  LM ; 
end 
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Appendix  B.  Reflective  Inverse  Diffusion  MATLAB 

Subfunctions 


function  CaptureMat rix  =  ReadSavelmage  (SBIG) 

%  Written  by  2Lt  Schafer,  adapted  from  code  written  by  Lt  Col  Mark 
%  Hoelscher  for  indirect  photography.  Updated  1  October,  2011. 

%%  Read  in  image  file  saved  by  SBIG  camera  and  return  as  a  matrix. 
fid=fopen (SBIG) ; 
fseek (fid, 2048, 'bof' ) ; 

Image=f read (fid,  [7  65,  510],  ' ushort ’ ) ; 

CaptureMat rix=Image ’  ; 
f close (fid) ; 

end 


function  Compact ImageMatrix  =  DeterminePixelValue (PV, Int , i, j , Compact ImageMatrix) 
%%  Return  the  compact  image  matrix  with  the  current  pixel  set  to  that  which 
%%  maximizes  the  target  intensity. 

[-i.  Index]  =  max  (Int); 

Compact ImageMatrix (i, j ) =PV (Index) ; 
end 


function  ImageMatrix=ExpandMat rix (StartMap, N) 

%%  Compute  number  of  pixels  of  the  SLM  grouped  into  one  effective  pixel. 
El  =  512/N; 

%%  Expand  compact  matrix  into  image  matrix  for  SLM 
for  col  =  1:N 
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for  row 


1  :N 


ImageMatrix (row*El— (El— 1)  : row*El,  col*El— (El— 1)  :col*El)  =  St artMap ( row, col ) ; 

end 

end 


function  Int  =  CalculatelntensityProf ileBeamSteer (Posit ion, .. . 

CaptureMat rix, PixRadius) 

%%  Written  by  Jessica  Schafer,  2Lt,  AFIT/ENP.  Updated:  1  Nov,  2011. 

%%  Adaptation  of  CalculatelntensityProf ile . m .  Adaptation  is  included  in 
%%  Appendix  B  instead  because  it  has  the  main  formula  with  extra 
%%  conditions . 

Spot lntensity=0 ; 

%%  Input  distances  between  center  spot  and  corner  spots  to  get  corner 
%%  target  intensity  regions, 
if  Position  ==  1 

CircDistVert ical  =  235; 

CircDistHorizontal  =  362; 
elseif  Position  ==  2 

CircDistVert ical  =  —235; 

CircDistHorizontal  =  362; 
elseif  Position  ==  3 

CircDistVert ical  =  235; 

CircDistHorizontal  =  —362; 
elseif  Position  ==  4 

CircDistVert ical  =  —235; 

CircDistHorizontal  =  —362; 

end 

for  j  =  1:765 

for  i  =  1:510 

if  (  (i<510/2+CircDistVertical&& j<7 65/2+CircDistHorizontal )  .  .  . 
&& (sqrt ( ( 510/2+CircDistVertical— i) "2+ (765/2+. . . 
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CircDist Horizontal— j ) "2 )  <PixRadius )  )  |  |  .  .  . 

( ( i>510/2+CircDistVert ical& & j<7 65/2+CircDistHorizontal ) . . . 

&& (sqrt ( ( 51 0/2+CircDistVert ical— (i— 1) ) "2+ 

( 7 65/2+CircDistHorizontal— j ) "2 )  <PixRadius )  )  |  |  .  .  . 

( (i<510/2+CircDistVertical&& j>7 65/2+CircDistHorizontal ) . . . 

&  & (sqrt ( (510/ 2+CircDistVert ical— i ) ~ 2+ ( 7  65/2  + 
CircDistHorizontal— ( j— 1 ) ) "2 )  <PixRadius )  )  |  |  .  .  . 

( ( i>510/2+CircDistVert ical& & j>7 65/2+CircDistHorizontal ) && . . . 
(sqrt ( (51 0/2+CircDistVert ical— (i— 1) ) ~2+(765/2+. . . 
CircDistHorizontal— ( j—1 ) ) "2 )  <PixRadius ) ) 

Spot Intensity=Spot Intensity+ . 5*CaptureMat rix (i,  j )  ; 

end 

if  (  (i<510/2+CircDistVertical&& j<7 65/2+CircDistHorizontal )  .  .  . 

&& (sqrt ( ( 510/2+CircDistVertical— i) "2+ (765/2+. . . 
CircDistHorizontal— j ) "2)  < (3*PixRadius / 4 )  )  )  1 1 .  .  . 

( ( i>510/2+CircDistVert ical& & j<7 65/2+CircDistHorizontal )  .  .  . 

&&  (sqrt  (  (51 0/2+CircDistVert  ical—  (i  —  1)  )  "2+(765/2+.  .  . 
CircDistHorizontal— j ) "2 )  < (3*PixRadius / 4  )  )  )  1 1 .  .  . 

( (i<510/2+CircDistVertical&& j>7 65/2+CircDistHorizontal ) . . . 

ScSc  (sqrt  (  ( 51 0/2+CircDistVert  ical— i )  "2+  (765/2  + 
CircDistHorizontal— (j  —  1) ) "2) < (3*PixRadius / 4  )  )  )  1 1 .  .  . 

( ( i>510/2+CircDistVert ical& & j>7 65/2+CircDistHorizontal ) . 
(sqrt ( (510/2+CircDistVertical— (i— 1) ) ~2+. . . 

( 7 65/2+CircDistHorizontal— ( j— 1) ) "2) < (3*PixRadius/4 ) ) ) 

Spot Intensity=Spot Intensity+ . 7  *CaptureMat rix (i , j ) ; 

end 

if  (  (i<510/2+CircDistVertical&& j<7 65/2+CircDistHorizontal )  .  .  . 

ScSc  (sqrt  (  ( 510/2+CircDistVert ical— i)  "2+  (765/2+.  .  . 
CircDistHorizontal— j ) "2 )  < (2*PixRadius / 4  )  )  )  1 1 .  .  . 

( ( i>510/2+CircDistVert ical& & j<7 65/2+CircDistHorizontal ) . . . 
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&  & ( sqrt ( (510/ 2+CircDistVert ical— (i— 1 ) ) "2  + (765/2+.  .  . 

CircDist Horizontal— j ) "2 )  < (2*PixRadius / 4 )  )  )  1 1  .  .  . 

( (i<510/2+CircDistVertical&& j>7 65/2+CircDistHorizontal ) . . . 

&&  (sqrt  (  (510/2+CircDistVertical— i)  ~2+(765/2+.  .  . 
CircDistHorizontal— ( j  —  1) ) "2) < (2*PixRadius / 4  )  )  )  1 1 .  .  . 

( ( i>510/2+CircDistVert ical& & j>7 65/2+CircDistHorizontal ) &&  .  .  . 
(sqrt ( (510 /2+CircDistVert ical— (i— 1) ) "2+(765/2+. . . 
CircDistHorizontal— ( j— 1) ) "2) < (2*PixRadius/ 4 ) ) ) 

Spot Intensity=Spot Intensity+ . 9*CaptureMat rix (i, j ) ; 

end 

if  (  (i<510/2+CircDistVertical&& j<7 65/2+CircDistHorizontal )  .  .  . 

ScSc  (sqrt  (  (510/2+CircDistVertical— i)  "2+  (765/2+.  .  . 
CircDistHorizontal— j ) "2 )  < (PixRadius / 4  )  )  )  1 1 .  .  . 

( ( i>510/2+CircDistVert ical& & j<7 65/2+CircDistHorizontal ) . . . 

Sc  Sc  (sqrt  (  (510/  2+CircDistVert  ical—  (i  —  1 )  )  "2  + (765/2+.  .  . 
CircDistHorizontal— j ) "2 )  < (PixRadius / 4  )  )  )  1 1 .  .  . 

( (i<510/2+CircDistVertical&& j>7 65/2+CircDistHorizontal ) . . . 

Sc  Sc  (sqrt  (  (510/2+CircDistVertical— i)  "2+(765/2+.  .  . 
CircDistHorizontal— ( j  —  1) ) "2)  < (PixRadius / 4 )  )  )  1 1 .  .  . 

( ( i>510/2+CircDistVert ical& & j>7 65/2+CircDistHorizontal ) . . . 

Sc  Sc  (sqrt  (  (510/  2+CircDistVert  ical—  (i  —  1 )  )  "2  + (765/2+.  .  . 
CircDistHorizontal— ( j  —  1) ) "2)  < (PixRadius/ 4 ) ) ) 

Spot Intensity=Spot Intensity+CaptureMatrix (i,  j)  ; 

end 


end 

end 

Int=Spot Intensity; 
end 


function  RecordAvgStart Intensities  (NMat ,  LUTMap, Folder ) 
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%%  Written  by  Jessica  Schafer,  2Lt,  AFIT/ENP.  Updated  1  Dec,  2011. 

%%  Determine  starting  average  speckle  intensity  with  random  phase  screens. 

ImageFrame  =  0; 

for  Iteration  =  1 : length (NMat ) 

N  =  NMat ( Iteration)  ; 
n  =  num2str (N,  ' % . 2d ’ ) ; 
for  Start  =  1:25 

S  =  num2str (Start )  ; 

StartMap  =  round (255*rand (N) ) ; 

ImageMat rix=ExpandMat rix (StartMap, N) ; 

BNS_LoadImageFrame  (ImageFrame, ImageMatrix, LUTMap) ; 

BNS_SendImageFrameToSLM (ImageFrame) ; 

ImageFrame=mod (ImageFrame  +  1,  26); 

Capturelmage; 

Init ialCapture  =  ReadSavelmage (' Capture . sbig ')  ; 

dlmwrite ( st rcat ( 1 C : \BNSMat lab SDK\ Schafer  SLM  Control \ ’, Folder, ' \ Initial ? , n, S, ’ .1 
InitialCapture,  ' delimiter '  ,  ’  \t ' ) ; 

imagesc (InitialCapture) ; colormap ( ' hot ’ ) ;  title (strcat ( ' Initial '  ,  n,  S) ) ; 

end 

end 

end 
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Appendix  C.  Surface  Profiles  of  Reflectors 


0.2  0.4  0.6  0.8  1 


Scan  Distance  (mm) 


Figure  33.  Surface  profiles  for  each  of  the  five  samples  measured  with  a  Tencor  Alpha 
Step  1000  profilometer.  The  RMS  surface  roughness  Rrms  is  given  for  each  sample. 
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Appendix  D.  Relative  Intensity  Plots  and  Images 


The  top  half  of  each  data  subset  displays  the  speckle  pattern  from  a  portion  of  the 
detector  center  for  phase  screens  of  N  =  4,  IV  =  8,  N  =  16,  and  N  =  32.  The  pattern 
is  normalized  to  the  average  starting  intensities  and  plotted  on  a  logarithmic  scale. 
The  bottom  half  shows  the  center  row  of  intensities  from  the  top  half  associated  with 
each  phase  screen.  Each  of  these  corresponds  to  the  final  iteration-not  necessarily 
the  most  optimal  iteration-for  its  respective  phase  screen  dimensions. 


Figure  34.  Relative  intensity  data  for  brushed  aluminum. 
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Figure  35.  Relative  intensity  data  for  sandblasted  aluminum. 
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Figure  36.  Relative  intensity  data  for  infragold. 
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Figure  37.  Relative  intensity  data  for  white  paint  on  reflection. 
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Figure  38.  Relative  intensity  data  for  graphite. 
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Figure  39.  Relative  intensity  data  for  spectralon. 
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